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ABSTRACT 


A  general  treatment  of  non-harmonic  vibrational  dynamics 
in  an  isolated  hydrogen  bond  is  presented,  based  on  an 
adiabatic  mode  separation  of  hydrogenic  motions  from  those 
of  the  heavier  particles. 

The  validity  of  the  Stepanov  (adiabatic)  approximation 
is  examined  by  including  correct  non-adiabatic  couplings  and 
solving  the  resultant  coupled  eigenvalue  problem  for  a  model 
system  (HF2-)  in  which  the  couplings  are  expected  to  have  a 
maximum  intrinsic  strength.  Convergence  is  demonstrated  by 
including  up  to  three  coupled  protonic  states.  Comparison 
with  exact  results  thus  obtained  shows  that  the  adiabatic 
approximation  gives  transition  frequencies  within  1%,  and 
infrared  relative  intensities  to  about  10%  -  in  cases  where 
the  coupled  protonic  levels  are  not  degenerate.  Cases  of 
isolated  strong  coupling,  arising  from  accidental  or  systematic 
degeneracies  of  proton  levels,  can  be  accommodated  in  the 
theory  in  a  direct  way,  but  do  not  arise  for  HF2-. 

Calculations  on  the  HF2-  system  are  based  entirely  on 
the  ab  initio  potential  energy  surface  of  Alml^f  (Ref.  24), 
but  employ  a  model  form  with  correct  dissociation  properties. 
The  close  agreement  of  computed  and  observed  transition 
frequencies,  and  IR  relative  intensities,  leads  to  a  convinc¬ 
ing  assignment  (in  the  adiabatic  mode  description)  for 
nearly  all  transitions  in  the  bifluoride  spectrum  (e.g.  in 
KHF2 ( s ) ) ,  on  the  basis  of  an  isolated  ion  model. 
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CHAPTER  ONE 


INTRODUCTION 


A.  Purpose. 

The  hydrogen  bond  has  been  the  subject  of  much 

experimental  and  theoretical  investigation.  This  has 

been  summarized  in  three  volumes  edited  by  Schuster, 

1  2  —  6 

Zundel,  and  Sandorfy,  in  several  other  texts,  and  in 

7 

many  review  articles  (surveyed  to  1977  by  Kollman  ) . 

The  presence  of  a  hydrogen  bond  (A-H..B)  is  strongly 
indicated  by  the  behaviour,  in  the  infrared  (IR)  spectrum 
of  the  band  associated  with  the  stretching  mode  v (A-H) . 
Upon  hydrogen  bonding,  this  band  is  lowered  in  frequency 

o 

and  greatly  broadened,  with  many  well  defined  peaks. 

In  addition,  the  frequencies  (^2'  v2 '  ^  Reformation 

modes  are  progressively  increased  with  increasing  hydrogen 
bond  (H-bond)  strength,  and  the  bands  are  also  generally 

o 

broad  and  structured.  It  is  now  widely  accepted  that 
this  broadening  and  structure  is  largely  due  to  strong 
anharmonic  coupling  between  the  A..B  symmetric  stretching 
mode  (v,)  ,  the  asymmetric  A-H  vibration  (v^)  ,  and  the 
deformation  modes  (v ^ t  v2 ' ^ *  In  Ueneral,  the  potential 
is  a  non-separable  function  of  the  four  vibrational 
coordinates,  and  cannot  be  made  separable  by  any  linear 
transformation . 
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2. 


However,  since  the  frequency  of  the  heavy  particle 
vibration  (v^)  is  invariably  much  less  than  the  frequencies 
associated  with  the  light  protonic  motions,  an  adiabatic 
mode  separation  was  suggested  by  Stepanov.  In  the 
Stepanov  approximation,  as  in  the  Born-Oppenheimer^ ' ^ 
separation  of  nuclear  from  electronic  motions,  a  set  of 
"protonic"  eigenstates  is  defined  for  each  configuration 
of  the  heavy  nuclei  (A...B).  These  states  depend  para¬ 
metrically  upon  the  "heavy  particle"  coordinate,  and 
each  protonic  eigenvalue  defines  an  effective  potential 
curve  for  the  heavy  particle  vibration  (v^) .  The  curves 
may  be  quite  different  for  each  protonic  state,  so  that 

much  of  the  breadth  and  complexity  of  the  spectrum  may 

12  13 

be  attributed  to  Franck-Condon  '  progressions  in  the 

g 

v-^  mode,  as  was  pointed  out  by  Sheppard. 

This  thesis  presents  a  general  theoretical  treatment 
of  the  vibrational  dynamics  of  an  isolated  H-bond,  using 
adiabatic  protonic  states  as  the  basis  for  a  representa¬ 
tion  of  the  exact  vibrational  wavefunctions .  In  this 
representation  (discussed  in  Chapter  Two)  there  results 
a  set  of  coupled  second  order  linear  differential  equations 
in  the  heavy  particle  coordinate,  the  solutions  to  which 
define  the  total  vibrational  eigenstates  of  the  system. 

The  Stepanov  approximation  is  a  "zero  order"  description 
in  this  representation:  the  non-adiabatic  coupling 
between  the  adiabatic  states  is  taken  to  be  negligible. 


1  I 

■ 
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3. 


It  is  necessary  to  distinguish  between  non-adiabatic 
coupling  and  the  "anharmonic  coupling"  which  arises 
because  the  H-bond  potential  function  is  non-separable . 

The  latter  is  not  negligible,  and  produces  a  strong 
dependence  in  the  protonic  potential  curves  on  the  heavy 
particle  coordinate,  and  a  reciprocally  strong  dependence 
in  the  dynamics  on  the  protonic  quantum  state.  These 
effects  may  be  accurately  described  in  the  Stepanov 
approximation  if  the  non-adiabatic  couplings  are  negligible. 

This  work  examines  the  effects  of  non-adiabatic 
couplings  on  the  validity  of  the  Stepanov  approximation. 
Calculations  have  been  performed  for  a  physically  realistic 
model  system  in  which  these  couplings  are  intrinsically 
as  strong  as  they  are  ever  likely  to  be  in  a  hydrogen 
bond.  The  conclusion  is  that  for  most  purposes  the  non- 
adiabatic  couplings  are  indeed  intrinsically  negligible. 

The  emphasis  on  intrinsic  strength  is  made  for  the 
following  reason:  in  cases  where  the  adiabatic  states 
are  degenerate,  or  nearly  so,  even  very  weak  non-adiabatic 
couplings  may  cause  the  states  to  be  strongly  mixed.  The 
corresponding  strong  vibronic  coupling  effects  are  well 
known  in  the  Born-Oppenheimer  separation  (e . g .  the  Jahn- 
Teller  effect).  In  hydrogen  bonds,  the  accidental 
degeneracy,  or  near  degeneracy,  of  the  protonic  levels 
'2v2'  and  *v3*  would  be  such  a  case;  the  strong  coupling 
("Fermi  resonance"14)  arising  in  this  situation  has  been 


. 


. 


4. 


.  .  15 

discussed  by  Witkowski  and  Wojcik.  In  this  thesis  the 
effects  of  non-adiabatic  couplings  are  examined  for  cases 
where  the  coupled  protonic  levels  are  not  degenerate. 

The  non-adiabatic  couplings  have  a  coordinate  system 

dependence.  In  the  usual  system  appropriate  to  a  molecular 

(H-bonded)  situation,  they  have  an  undesirable  property: 

they  do  not  vanish  as  the  H-bond  dissociates.  This  defect 

has  been  corrected  by  applying  the  non-linear  coordinate 

1 6 

transformation  technique  of  Thorson  and  Delos  (Chapter 
Two,  section  D3) .  A  degree  of  approximation  is  introduced, 
but  the  dominant  effects  of  extending  a  calculation 
beyond  the  Stepanov  description  are  retained. 

The  system  used  to  examine  non-adiabatic  coupling 
effects  is  a  realistic  model  for  the  (v^,  v^)  stretching 
dynamics  in  the  strongly  hydrogen-bonded  bifluoride  ion 
(HF2-).  Calculations  were  carried  out  using  a  model 
potential  surface  based  entirely  on  ab  initio  calculations. 
This  model  is  described  in  Chapter  Three. 

It  is  also  necessary  to  establish  techniques  for  the 
accurate  numerical  solutions  of  (a)  the  Schroedinger 
equation  describing  essentially  proton  motion  at  a  fixed 
heavy  nuclear  configuration  (yielding  "vibrational  Born- 
Oppenheimer"  states,  Chapter  Four),  and  (b)  the  set  of 
coupled  equations  (Chapter  Five) .  Techniques  which  are 
applicable  to  a  linear  symmetric  system  (A-H...A)  are 
described.  This  is  not  a  restriction  imposed  by  the 


8 
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5. 


general  theory,  but  it  does  allow  a  great  reduction  in 
computation  time. 

For  HF2~ ,  the  calculated  IR  and  Raman  transition 
frequencies,  and  IR  relative  intensities,  allow  assignments 
of  lines  in  the  experimental  spectrum  of  crystalline  KHF2 
to  quantum  states  defined  within  the  Stepanov  description 
(Chapter  Six) .  This  crystal  spectrum  can  be  well  under¬ 
stood  using  the  isolated  ion  model. 

The  theoretical  calculation  shows,  quite  independently 
of  the  experimental  agreement,  that  with  an  adiabatic 
separation  (Stepanov  description)  the  useful  concept  of 
mode  excitations  may  be  retained,  while  accurately  includ¬ 
ing  the  effects  of  a  non-separable ,  strongly  anharmonic 
interaction.  However,  many  of  the  predictions  of  a 
simple  harmonic  mode  description  must  be  abandoned: 
overtone  frequencies  are  not  integer  multiples  of  funda¬ 
mentals,  the  frequency  of  the  heavy  particle  mode  (v^) 
depends  strongly  on  the  protonic  state,  and  intensities 
may  differ  greatly  from  those  predicted  using  a  harmonic 
basis;  yet  the  common  idea  is  that  each  level  can  be 
assigned  distinct  quantum  numbers  for  each  of  a  collection 
of  separate  modes.  A  simple,  yet  general  picture  of  the 
vibrational  dynamics  of  an  isolated  hydrogen  bond  thus 


emerges . 


' 
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B.  Previous  Vibrational  Calculations. 

1.  General  Formulation. 

17 

Janoschek  has  reviewed  the  methods  which  have  been 
used  in  vibrational  calculations  to  date.  The  usual 
procedure  is: 

(a)  centre  of  mass  motion  is  extracted  and  ignored; 

(b)  the  Born-Oppenheimer  separation  of  nuclear  and 
electronic  motions  is  assumed  to  be  valid; 

(c)  vibration-rotation  interactions  are  ignored; 

(d)  the  equilibrium  configuration  is  taken  to  be 
linear . 

At  this  point  the  internal  motion  of  the  three  nuclei 
(A-H..B)  may  be  expressed  in  terms  of  four  vibrational 
coordinates  which  correspond,  in  the  harmonic  approxima¬ 
tion,  to  "proton  stretching"  and  "A...B  stretching"  modes, 
and  the  degenerate  "bending  modes".  A  further  approxima¬ 
tion  is  then  introduced: 

(e)  bending  and  stretching  modes  are  assumed  to  be 
separable . 

1 8 

Coulson  and  Robertson  have  pointed  out  that,  since 
accurate  experimental  frequencies  for  the  bending  modes 
are  not  generally  available,  it  is  of  interest  to  see  what 
progress  can  be  made  by  ignoring  them. 

In  the  absence  of  an  accurate  four-dimensional 
potential  function,  some  approximate  separation  of  bending 
and  stretching  is  clearly  necessary. 


In  Chapter  Two  it  is  assumed  that  the  bending  and 
stretching  motions  of  the  proton  are  approximately  separ¬ 
able  at  any  fixed  heavy  nuclear  separation  (R) .  The 
separate  potentials  are  then  both  taken  to  be  parametrically 
dependent  on  R.  It  follows  that  the  adiabatic  curves  for 
heavy  nuclear  motion  are  dependent  upon  both  bending  and 
stretching  "vibrational  Born-Oppenheimer "  eigenstates  of 

the  proton.  Any  Fermi  resonances  between  bending  and 

15 

stretching  modes  of  the  proton  can  be  treated  by 
introducing  their  interaction  as  a  perturbation  within 

the  manifold  of  near-degenerate  proton  states. 

19 

There  is  evidence  that  the  R-dependence  of  the 

potential  function  for  is  much  weaker  than  that  for 
Vy  which  implies  that  an  approximate  treatment  ignoring 
bending  interactions  may  yield  useful  results.  This  is 
the  approach  taken  here  for  the  specific  study  of  the 
bifluoride  ion. 

The  main  focus  of  attention,  in  past  calculations 
and  in  this  thesis,  is  then  upon  the  two  stretching 
coordinates,  and  the  non-separable  two-dimensional  potential 
surface  which  describes  their  interaction. 

If  a  pair  (q^,  q^)  of  generalized  Jacobi  coordinates 
(discussed  in  detail  in  Chapter  Two)  is  used  (correspond¬ 
ing  to  the  stretching  modes  (v-^,  v^)  /  with  associated 
reduced  masses  (m1 ,  m3) ) ,  the  vibrational  Schroedinger 
equation  contains  no  cross-terms  in  the  kinetic  energy 


V1 


. 


8. 


operator : 

{-(h2/2m1)3  2/3q2  -  (fl2/2m3 )3  2/3  q2  +  V  (q3  ^  -E } ^  (q3  =  0 

(1.1) 

The  solution  of  (1.1),  with  V(q^fq^)  from  ab  initio 

calculation,  or  from  empirical  data,  or  both,  has  been 

17 

the  main  subject  of  recent  vibrational  calculations. 

The  methods  which  have  been  used  fall  into  two  classes: 

(a)  "direct"  solutions  of  the  two-dimensional  equa¬ 
tion  by  expansion  in  products  of  harmonic 
oscillator  functions; 

(b)  adiabatic  representations. 

A  brief  discussion  of  these  methods  follows. 

2.  Direct  Two-Dimensional  Solutions. 

An  eigenfunction  ) /  which  satisfies  (1.1) 

under  appropriate  boundary  conditions,  may  be  expressed 
as  a  linear  combination  of  binary  products  of  harmonic 
oscillator  (ho)  functions  in  q^  and  q^: 

^n(q3'ql}  =  £  Cnii *Xi (q3} *Xi (ql}  ;  (1*2) 

if  j 

where 

Xi  (qk)  =  Hj_  (nkqk)  *exp  (-0.5n^*qk2)  ,  (k  =  1  or  3)  (1.3) 

and  {H^}  are  normalized  Hermite  polynomials. 

The  coefficients  {Cn^j}  are  obtained,  for  a  given 
basis  size  N,  by  application  of  the  variational  principle. 


9. 


Diagonalization  of  the  resulting  NxN  matrix  yields  upper 

bounds  on  the  first  N  eigenvalues.  The  exponents  (n^f  n^) 

may  initially  be  fixed  using  either  the  transition 

20  21 

frequencies  and  (as  Ibers,  and  Singh  and  Wood  ), 

or  the  curvature  at  the  minimum  of  the  potential  surface 

22  23 

(as  Janoschek  et  al .  '  ) ,  or  they  may  be  treated  as 

24 

further  variational  parameters  (as  Alml0f  ).  In  principle, 
the  values  for  {ri^}  are  not  critical,  because  the  functions 
{X^  (q^.)  ;  i  =  1,  2,  .  ..»}  are  complete  in  q^  space,  but 
a  judicious  choice  will  speed  the  convergence  of  (1.2). 

The  problems  of  convergence  and  interpretation  are 
the  main  difficulties  with  this  approach  to  the  solution 
of  (1.1). 

Without  exception,  those  who  have  applied  this 
"matrix  expansion"  method  have  used  a  polynomial  form 
for  the  representation  of  the  potential  surface: 

V(q3,qi)  =  l  akJ,q*q*  U-4) 

k ,  £ 

The  coefficients  {a  }  are  then  obtained  either  from 

20  21 

experimental  data,  '  or  by  least-squares  fit  to  ab 

22-24 

initio  data  points.  The  number  of  terms  used  m 

20  .  . 

(1.4)  varies  greatly:  the  form  of  Ibers  is  limited  to 

four  (with  anharmonic  coupling  arising  only  through  a 
2  24 

term  in  q^q^) ,  whereas  Alml^f  included  eight  terms. 

The  polynomial  representation  for  the  potential 
function  is  used  with  good  reason:  matrix  elements  {H^}, 
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arising  in  the  secular  equations,  are  readily  calculated. 

Unfortunately,  it  has  been  found  necessary  to  use  large 

basis  sets  in  order  to  ensure  convergence  for  the  lower 

21 

eigenvalues:  Singh  and  Wood  required  N  =  144,  and 


24 

Alml^f  had  N  =  160. 

22  23 

Janoschek  et  at.  '  have  used  non-orthogonal  basis 
sets  for  double-minimum  surfaces  (with  the  ho  functions 
localized  at  the  minima:  a  two  centre  expansion),  in 
order  to  provide  more  rapid  convergence  (N  =  32) .  This 
introduces  overlap  integrals  {S^ ^ }  into  the  secular 
equations,  and  consequently  the  possibility  of  over¬ 
completeness  for  the  higher  eigenstates  (as  N  increases, 
or  as  the  minima  coalesce,  the  matrix  becomes  singular) 

The  question  of  convergence  to  an  accurate  eigen¬ 
value  is  particularly  disturbing  with  a  polynomial 
representation  for  the  potential  surface.  The  surface 
given  by  (1.4)  may  be  accurate  in  the  region  of  the 
equilibrium  configuration  (where  ab  initio  or  empirical 
data  have  been  used  to  define  the  coefficients) ,  but  as 
the  hydrogen  bond  dissociates  the  representation  inevit¬ 
ably  becomes  very  poor. 

In  addition,  having  achieved  convergence  with  many 
ho  products,  it  may  become  difficult  to  interpret  the 
associated  eigenfunctions  in  a  conceptually  useful  manner 


Previous  vibrational  calculations  for  the  bifluoride 
ion20'24  have  used  this  approach,  and  will  be  further 
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discussed  in  Chapter  Three.  Other  systems  which  have 
been  investigated  using  this  two-dimensional  treatment 
are  H5C>2+,  22'25  and  H3C12±.  23 

21 

The  work  of  Singh  and  Wood  is  of  particular 
relevance  to  this  thesis,  as  they  have  made  a  quantitative 
estimation  of  the  validity  of  the  Stepanov  approximation 
for  the  lower  lying  levels  of  a  model  H-bonded  system. 
Their  four  term  polynomial  representation  of  the  potential 
function  (Ibers  form  )  was  chosen  to  approximate  the 
symmetric  double  minimum  function  appropriate  to  0-H..0 
systems.  Coefficients  a^Q,  aQ^  and  aQ^  (equation  (1.4)) 
were  fixed  by  requiring  175  cm  ^ 3000  cm  ^ , 

and  a  barrier  height  5000  cm  \  while  the  anharmonic 
coupling  parameter  a^2  was  related  to  an  equilibrium 
configuration.  Their  results  show  very  close  agreement 
(within  1%  for  the  lower  levels)  between  eigenvalues 
calculated  by  the  two-dimensional  matrix  expansion  method, 
and  the  corresponding  Stepanov  approximation  values 
(which  they  obtained  by  first  solving  the  one-dimensional 
Schroedinger  equation  (ho  basis  in  q^)  for  m^ ,  at  fixed 
q^  values,  yielding  protonic  eigenvalue  curves,  and  then 
solving  the  one-dimensional  equations  (ho  basis  in  q^) 
using  these  curves  as  potential  surfaces  for  m^) . 

With  this  model  they  have  convincingly  shown  that 
the  Stepanov  approximation  is  very  good  for  H-bonds  of 
weak  to  medium  strength,  involving  for  example,  oxygen 


. 


and  nitrogen,  but  they  have  not  quantitatively  investigated 
the  energy  levels  for  a  strong  bond  such  as  the  bifluoride 
ion . 

2  g 

In  a  later  paper  Singh  and  Wood  have  studied  the 

Stepanov  approximation  for  the  bifluoride  ion,  but  only 

in  relation  to  the  contraction  (AR)  in  equilibrium  F-F 

2  0  2  8  —  3  0 

distance  upon  deuteration.  '  Again  using  a  four 

term  polynomial  representation  (with  the  coefficients  of 

27 

Ibers  and  Delaplane,  empirically  estimated  for  HF2~) 
for  the  potential  function,  they  found  that  a  Stepanov 
calculation  accurately  reproduces  the  isotope  contraction, 
and  have  suggested  that  the  lower  energy  states  are 
therefore  well  approximated. 

3.  Adiabatic  Representations. 

The  Hamiltonian  of  (1.1)  may  be  written: 

H(q3,q1)  =  (p^/2m1)  +  Mq^q^  ,  (1.5) 

where 

h(q3;q1)  =  (p3/2m3)  +  V(q3;q1)  , 

and  Pj  is  the  momentum  operator,  -ih-9/8  q ^  (j  =  1,3). 

Recall  that  q1  and  q3  are  a  pair  of  generalized 
Jacobi  coordinates  with  associated  relatively  heavy  and 
light  reduced  masses  (n^  and  m3)  respectively.  There 
are  three  possible  choices  for  this  pair  (described  in 


(1.6) 
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Chapter  Two) ,  each  of  which  is  appropriate  for  a  particular 
configuration  (A-H....B,  A-H-B,  A....H-B).  There  is  an 
important  significance  to  this  choice,  which  will  be 
discussed  in  detail  in  Chapter  Two,  in  connection  with 
the  asymptotic  behaviour  of  the  non-adiabatic  coupling 
terms.  For  any  choice,  however,  the  adiabatic  representa¬ 
tion  employs  eigenfunctions  of  the  "protonic"  Hamiltonian 
hCq^jq-^)  (calculated  at  any  fixed  value  of  q^)  as  basis 
states  for  an  expansion  of  the  total  wavef unction : 

N 

’Vq3'ql)  =  l  fj£<ql><^£<q3'ql)  '  (1-7) 

X/ 

with  {(J)^ }  defined  by: 

h  (q3;q1)  (q3;q1)  =  u£  (SiHjj,  (q3;qx)  d*8) 

It  is  shown  in  Chapter  Two  that  substitution  of  (1.7) 
into  (1.1)  results  in  an  NxN  set  of  coupled  equations 
for  the  functions  {f in  q±  space.  The  coupling  arises 
from  terms: 

pk*(ql»  5  /+”  *k  Pl*ldq3  '  (1‘9) 

—  00 

The  adiabatic  (Stepanov)  approximation  neglects  these 
terms  (i.e.,  assumes  the  protonic  states  change  "slowly" 
with  respect  to  q1) .  The  equations  which  result  in  the 
adiabatic  approximation  are  then  very  simple: 

lPl/2ml  +  Vql)]fj£(ql>  =  Ej£fj£(ql>  ‘ 


(1.10) 


. 
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Equation  (1.10)  describes  the  motion  of  a  mass  m^ 
on  a  proton-state-dependent  potential  curve  (£  gives 
the  state  of  excitation  of  the  proton) .  For  each  protonic 
curve  there  are  associated  eigenvalues  {E . 0 ;  j  =  0,  1, 

2,  ...},  and  so  the  adiabatic  separation  has  retained  the 
concept  of  a  mode  description :  total  eigenstates  may 
be  labelled  by  a  double  index:  |j,£>.  Changes  in  £  (for 
fixed  j)  refer  to  excitation  of  the  "A-H  stretching  mode" 

(v3) ,  and  changes  in  j  (for  fixed  £)  refer  to  excitation 
of  the  "A..B  mode"  (v^)  .  Of  course,  these  are  not 
excitations  of  any  "pure  mode"  in  the  sense  of  a  harmonic 
separation,  but  it  is  a  useful  concept.  Because  of  a 
shift  in  the  minima  and  form  of  the  excited  state  curves 
{U^},  a  Franck-Condon  effect  is  expected:  the  combination 
bands  "mv^"  i  "nv^"  will  have  significant  intensities. 

The  adiabatic  description  has  been  extensively 

employed  by  Marechal,  Witkowski,  and  numerous  co- 
15  31-45 

workers.  '  As  they  were  concerned  with  asymmetric 

H-bonds  of  weak  to  medium  strength,  they  have  treated  the 
potential  function  in  the  following  empirical  manner: 
the  proton  motion  is  assumed  to  be  harmonic,  but  with 
force  constant  and  equilibrium  position  parametrically 
dependent  upon  the  A...B  separation.  They  write: 

V(q,Q)  =  l/2-mw2(Q)  [q-qQ(Q)]2  +  l/2-Mfi2(Q-Q0)2  (1.11) 

(where  q  is  to  be  identified  with  q3,  Q  with  q1#  and  (m,M) 
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with  (m^m^);  w  and  Q  are  protonic  and  heavy  particle 
frequencies  respectively) . 

By  expanding  w(Q)  and  (Q)  about  the  equilibrium 
value  Qq,  in  powers  of  Q,  and  retaining  only  linear 
terms,  they  have  related  the  principal  anharmonic  coupling 
to  the  parameter  (dw/dQ)^  .  By  a  quantitative  examination 

Qo 

46 

of  the  potential  surfaces  of  Reid  (which  are  based  on 

the  empirical  potential  function  of  Lippincott  and 
4  7 

Schroeder,  for  0-H. . .0  systems)  they  concluded  that 
(dqQ/dQ)Qo  and  higher  terms  are  of  much  less  importance. 

In  their  model  the  ground  state  protonic  curve  UQ (Q) 
is  harmonic: 

Uq(Q)  -I  const*  +  1/2  •  (Q-Qq)  ^  ,  (1.12) 

whereas  the  curve  for  the  first  excited  state  of  the 
proton  is  vertically  displaced  and  linearly  distorted: 

U1(Q)  ~  Uq(Q)  +  r  +  bQ  ,  (1.13) 

15 

where  r  is  the  "vertical  excitation  energy"  (essentially 
v 3 ) ,  and  b  is  related  to  (dw/dQ) qq .  In  practice  b  is 
altered  to  give  a  good  spectral  fit  (it  is  an  empirical 
parameter  in  their  theory) . 

3  3 

Within  this  model,  they  have  estimated  the  effect 
of  non-adiabatic  coupling  as  a  perturbation  of  harmonic 
functions  in  Q  space.  They  find  the  perturbation  energy 
to  be  about  1/300  of  the  proton  quantum  energy,  hw(QQ), 
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and  consequently  discard  the  non-adiabatic  coupling  terms. 

Assuming,  then,  the  validity  of  the  Stepanov  approx¬ 
imation  within  an  isolated  H-bond,  and  introducing 
empirical  forms  for  the  surfaces  {U^},  their  work  then 
explores  couplings  between  identical  A-H..B  units,  (e.g., 
in  carboxylic  acid  dimers  )  which  have  degenerate 
Stepanov  eigenstates  in  zero  order. 

Similar  adiabatic/harmonic  models  have  been  used  to 
study  interactions  between  many  bonds  (e.g.,  in  the 
lattice),  and  have  been  reviewed  by  Hofacker,  Marechal, 
and  Ratner.^ 

A  similar  approach  is  also  taken  by  Sokolov  and 
49 

Savel'ev,  who  have  shown  that  the  main  characteristics 
of  the  H-bond  spectra  of  asymmetric  (A-H..B)  and  symmetric 
(A-H. .A)  systems  may  be  understood  within  the  adiabatic 
approximation.  They  have  also  estimated  the  adiabatic 
coupling  effects  as  a  perturbation  within  the  adiabatic/ 
harmonic  model,  and  have  concluded  that  even  for  the 
strongly  bonded  bifluoride  ion,  the  effects  are  small  for 
the  low  energy  states. 

50  .  . 

Coulson  and  Robertson  have  assumed  the  validity 

of  the  adiabatic  approximation  in  their  detailed  study 
of  MeOHX  (X  =  Cl,  F)  systems.  Their  potential  function 
model  is  very  simple,  consisting  of  the  sum  of  either  a 
quadratic  or  a  Morse  function  in  the  X-H  bond  length 
(they  consider  both  cases) ,  and  a  Morse  function  in  the 


. 

if 


0-H  distance.  Even  with  such  a  simple  function,  for 

which  the  anharmonicity  of  the  X-H  stretching  mode  is 

unaffected  by  H-bond  formation,  the  potential  is  non- 

separable  in  an  appropriate  Jacobi  coordinate  system, 

and  gives  rise  to  substantial  combination  band  intensities 
51 

Robertson  extended  the  calculation  to  include  a 
potential  term  explicitly  dependent  upon  both  X-H  and 
0-H  distances  (effectively  allowing  the  potential  for 
the  X-H  mode  to  be  X-0  dependent) ,  but  he  still  calculated 
line  positions  and  relative  intensities  within  the 
adiabatic  approximation. 

52 

Finally,  Anderson  and  Lippmcott,  and  Jiang  and 

53-55  47 

Anderson  (using  modified  Lippmcott  and  Schroeder 

forms  for  the  potential)  have  numerically  calculated  the 

protonic  eigenvalue  curves  for  several  systems.  They 

estimated  "vertical  transition  frequencies"  between  these 

curves,  but  have  not  discussed  non-adiabatic  couplings. 

53 

The  functional  form  of  Jiang  and  Anderson  for  the 
potential  in  HF2~,  has  been  adapted  in  this  thesis  to 
give  a  reasonable  fit  to  the  ab  initio  data  points  of 
Alml^f,2^  and  will  be  discussed  in  detail  in  Chapter  Three 
Evidently  the  adiabatic  separation  of  heavy  nuclear 
and  protonic  motions  is  widely  accepted  in  the  literature, 
for  both  weak  and  strong  hydrogen  bonds.  However,  none 
of  these  studies  in  the  adiabatic  representation  have 
explicitly  investigated  the  form  of  the  non-adiabatic 


■ 
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couplings,  or  have  included  them  in  any  extension  of  a 
calculation  beyond  the  Stepanov  approximation. 


CHAPTER  TWO 


GENERAL  THEORY 

A.  Vibrational  Coordinates. 

1.  Centre  of  Mass  Motion. 

For  a  system  of  N  particles  of  total  mass  M,  upon 
which  no  external  forces  act,  the  motion  of  the  centre 
of  mass  (CM)  is  separable  from  the  relative  motions  of 
the  particles,  and  is  just  the  free  translation  of  a 
particle  of  mass  M.  In  the  following  treatment  this 
motion  has  been  extracted  and  ignored  (in  a  crystal 
environment  (e.g.,  KHF2)  this  degree  of  freedom  is 
associated  with  lattice  modes) .  The  internal  motion  of 
the  system  is  then  described  using  N-l  relative  vectors. 

2.  Jacobi  Relative  Coordinates. 

The  use  of  N  generalized  Jacobi  vectors,  one  of 

which  is  the  vector  position  of  the  CM,  leads  to  a  diagonal 

kinetic  energy  operator  (containing  no  cross  terms) .  A 

5  6 

detailed  account  is  given  by  Blokhintsev.  For  three 
particles  (nine  degrees  of  freedom) ,  a  Jacobi  system 
consists  of  3  external  coordinates  (defining  the  position 
of  CM) ,  and  6  internal  coordinates  (two  relative  vectors) . 
One  internal  vector  connects  any  chosen  pair  of  the 
particles,  and  the  other  gives  the  position  of  the  third 
particle  from  the  centre  of  mass  of  the  chosen  pair. 
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The  definition  of  a  Jacobi  system  does  not  specify 
any  preferred  pair,  and  the  three  different  sets  are 
each  appropriate  in  a  different  physical  configuration 
(section  D2  of  this  Chapter) . 

3.  Separation  of  Vibration  and  Rotation. 

Under  the  assumption  that  the  vibration  of  a  molecular 

system  does  not  significantly  alter  the  equilibrium  inertial 

axes  and  moments  (i.e.,  the  molecular  frame  is  essentially 

rigid)  the  approximate  separation  of  rotation  and  vibra- 

57 

tion  is  well  established.  If  the  H-bond  is  not  isolated 
the  rotational  degrees  of  freedom  may  appear  as  lattice 
modes  (e.g.,  in  the  bifluoride  salts),  or  as  "rocking 
modes"  (e.g.  in  the  gas  phase  ether-HCl  systems'^  ^) . 
Interactions  with  such  modes  will  not  be  considered  here. 

It  will  be  assumed  that  the  H-bond  has  a  linear  equilibrium 
configuration,  and  the  free  rotation  of  the  inertial  axis 
will  be  removed  and  ignored.  This  axis  then  defines  the 
reference  frame  for  the  four  (3N-5)  vibrational  coordinates. 

Using  Jacobi  coordinates  R  and  r  (Figure  2.1) 
appropriate  to  the  "molecular"  configuration,  the  vibra¬ 
tional  Schroedinger  equation  is: 

'  {-(h2/2y)9  2/3R2  -(h2/2m)^r2  +  V(r,R)  -E}ij;(r,R)  =  0  ;  (2.1) 

where  y  is  a  "heavy  particle"  reduced  mass: 

y  =  MaMb/(Ma  +  Mb> 


I 


(2.2) 


>  cp 


FIGURE  2.1 


MOLECULAR  VIBRATIONAL  COORDINATES 


CMN  =  CENTRE  OF  MASS  OF  A, B 
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and  m  is  the  "molecular"  proton  reduced  mass: 


m  =  Mr(Ma  +  Mb)/Mt 


(2.3) 


Ma,  Mg,  Mg  are  the  atomic  masses,  and  MT  is  their  sum. 
The  solution  of  (2.1),  with  appropriate  boundary 


conditions,  yields  vibrational  eigenvalues  {E ^ }  and 
corresponding  eigenfunctions  (r,R) }. 


4.  Harmonic  and  Adiabatic  Mode  Descriptions. 

One  may  strictly  refer  to  isolated  vibrational  modes 
only  if  the  potential  surface  V(r,R)  is  a  separable  func¬ 
tion  of  the  four  coordinates.  In  such  a  case  there  would 
be  two  stretching  modes  (symmetric  and  asymmetric) ,  and 
two  bending  modes  (degenerate  for  an  isolated  linear 
bond) .  These  modes  might  be  harmonic  in  a  first  approxima¬ 
tion  . 

The  H-bond  potential  is  certainly  non-separable ;  but 
if  an  adiabatic  separation  is  shown  to  be  possible,  then 
a  mode  description  may  nevertheless  be  retained  (as 
discussed  in  Chapter  One,  section  B3) . 


B.  Adiabatic  Representation. 

1.  Expansion  Basis. 

If  at  any  point  R^ ,  in  R  space,  there  exists  a  complete 
set  of  functions  '  {(f>£  (r;R±)  }  which  span  the  r  space,  then 
any  function  of  r,  defined  at  R^ ,  may  be  expanded  in  this 
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set.  With  the  assumption  of  continuity  in  R,  for  each 
member  of  {(f)^},  the  wavefunction  (r,R)  may  be  written: 

<|)j(r,R)  =  l  f  .  £  (R)  ‘  <t>£  (r ;  R)  .  (2.4) 

.2/ 

In  the  adiabatic  representation,  the  complete  set 
{(f)^}  is  chosen  to  be  the  set  of  eigenfunctions  of  the 
Hamiltonian  for  the  reduced  mass  m  (i.e.,  essentially 
protonic  vibrational  Born-Oppenheimer  states) : 

h  (r ;  R)  c})^  (r ;  R)  =  U£  (R)  (r ;  R)  ,  (2.5a) 

where 

h(r;R)  =  -(h2/2m)Vr2  +  V(r;R)  .  (2.5b) 

In  equations  (2.5),  the  semi-colon  indicates  that  R 
is  treated  as  a  fixed  parameter.  At  any  R,  the  functions 
are  eigenfunctions  of  an  Hermitian  operator  in  r, 
and  therefore  form  a  complete  set  in  r,  with  corresponding 
eigenvalues  {U^}. 

2.  Coupled  Equations. 

Using  the  expansion  (2.4),  with  {4>^}  defined  by  (2.5), 
the  Schroedinger  equation  (2.1)  gives: 

J{- (Ti2/2p )  32/3R2  +  [U£(R)-E.]  }f  .£(R)iJ>£(r;R)  =  0  .  (2.6) 

l 

For  any  operator  A,  the  bracket  notation  shall  be 


defined  by: 


an oi  tomtf  .1$  .  *  **  al 


- 


24. 


<k  |  A  |  £>  E  f  0*  (r  ;  R)  A(j)  ^  (r  ;  R)  dr 


(2.7) 


so  that  (2.6)  becomes: 

I{-(ft2/2,)[<k|a2/3R2|*>fjjl(R)  +  f  (R)  (2.8) 

+  2<k|3/3R|£>f’£  (R)]  +  [U^tRJ-Ejlf  .£(R)6kJ,}  =  0  . 

A  coupling  matrix  P(R)  is  now  defined,  with  elements: 
Pk£(R)  E  “ih<kl9/9RU>  /  (2.9) 


which  leads  to  the  identity: 

-Ti2<k  |  32/3R2  |  «■>  =  -ihP^  (R)  +  [P2  (R)  ]  k£  ;  (2.10) 


and  (2.8)  may  be  written: 


|{-h2/2p)[(i/ft)P-Jl  .fjt  -  d/ft2)  (P2)kjl-f  j4  +  fgt  -«k4 

+  (2i/h)Pk£-f^  ]  +  [U£-Ej  ]  f  ^  *  6k£}  =  0 


(2.11) 


This  takes  the  simple  matrix  form: 


{ (l/2y) [-ih-d/dR-I  +  P(R) ]2  +  [U  (R) -E..  I]  }f  (R)  =0  (2.12) 


f.  (R)  is  a  column  vector  with  components  {f..}, 

J  J 

I  is  the  unit  matrix,  and  U (R)  is  diagonal  with  elements 


{V- 


These  coupled  equations  were  originally  derived,  for 

6 1 

a  general  adiabatic  representation,  by  Born  and  Huang. 

In  principle,  an  exact  solution  to  the  vibrational 
eigenvalue  problem  may  be  obtained  by  solving  the  two  sets 
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of  equations,  (2.5)  and(2.12),  with  sufficient  terms  in 
the  expansion  (2.4). 

3.  Stepanov  Approximation. 

The  matrix  equation  (2.12),  resulting  from  the 
adiabatic  representation,  shows  that  negligible  coupling 
situations  (P  *  0)  lead  to  the  Born-Oppenheimer  separa¬ 
tion  of  fast  and  slow  relative  motions.  The  motion  of 
the  "slow"  heavy  reduced  mass  y  is  then  determined  by 
one  of  the  protonic  energy  curves  {U^ (R) }  through  the 
uncoupled  differential  equation: 

{- (ft2/2U)d2/dR2  +  (R)  -  Ej^Jf  (R)  =  0  (2.13) 

These  equations  (one  for  each  eigenstate  of  the 
proton:  l  =  0,  1,  2,  ...)  provide  the  basis  for  the 

Stepanov  approximation  (cf.  equation  (1.10)). 

In  this  case  the  total  eigenfunctions  have  a  single 
product  form: 

ijj  .  ^  (r  ,R)  =  f  j  ^  (R)  *  <J>£  (r ;  R)  (2.14) 

(one  term  in  expansion  (2.4)),  and  together  with  eigen¬ 
values  { E  •  0 } ,  may  be  indexed  by  to  indicate  the 

protonic  state  with  which  they  are  associated. 

A  description  invoking  excitation  of  the  "protonic 
modes"  v3  and  v2  (for  which  A  changes),  and  the  "A...B 
stretching  mode"  v1  (for  which  j  changes)  may  be  retained, 
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and  their  strong  interaction  will  be  indicated  by  shifts 
in  the  minima  and  "force  constants"  of  the  curves  {UQ}. 

The  intensities  of  "combination  bands"  (e.g.,  v3  ±  v1) 
may  be  significant,  as  predicted  by  the  Franck-Condon 
principle . 

The  validity  of  the  Stepanov  approximation  is  a  major 
concern  of  this  thesis. 

If  the  coupling  between  one  or  more  degenerate  or 
nearly  degenerate  levels  is  not  zero,  then  the  Stepanov 
approximation  breaks  down,  and  the  wavef unction  must  be 
represented  by  an  expansion  (2.4),  with  the  sum  including 
all  the  (near- ) degenerate  states.  An  important  case  is 
the  possible  near-degeneracy  of  the  stretching  funda¬ 
mental  (v3)  and  the  bending  overtone  (2v3) .  This  case 

15 

has  been  discussed  by  Witkowski  and  Wojcik. 

This  possible  degeneracy  represents  a  "special  case" 
in  the  Stepanov  approximation,  and  does  not  occur  in  the 
HF3-  system  considered  here  as  a  test  problem  for  the 
"normal"  validity  of  the  approximation. 

C.  Separation  of  Bending  and  Stretching  Potentials. 

1.  Approximations. 

In  attempting  a  solution  of  (2.5): 

{- (fi2/2m)  Vr2  +  V(r;R)  }  i}>  £  (r;R)  =  (R)  <t>^  (r ;  R)  ,  (2.5) 

the  following  approximate  separation  of  the  potential 


■ 


■ 
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will  be  made: 

V  (r ;  R)  ~  Vz(z;R)  +VKy(x,y;R)  ,  (2.15) 

where  the  z  coordinate  of  r  is  defined  to  be  collinear 
with  R. 

For  an  isolated  linear  H— bond,  the  potential  energy 
will  not  be  dependent  on  the  cylindrical  polar  angle  <J> 
(cylindrical  symmetry),  and  (2.15)  may  be  written: 

V (r ; R)  V  (z;R)  +V  (p;R)  (2.16a) 

^  p 

2  2 

where  p  =  x  +  y  . 

It  may  be  further  assumed  that  Vp  is  approximately 
harmonic,  but  with  an  R-dependent  force  constant,  giving: 

V (r ; R)  ~V  (z;R)  +  1/2-K  (R)p2  .  (2.16b) 

z  p 

With  (2.16b)  the  Schroedinger  equation  is  separable 

in  cylindrical  polar  coordinates.  The  details  are  given 

6  2 

by  Pauling  and  Wilson. 

The  bending  motion  has  been  reduced  to  that  of  a 
two-dimensional  isotropic  oscillator  in  the  x,y  (p,<|>) 
plane . 

[If  it  is  observed  (V2  f  ^2 ' )  that  environmental 

effects  lift  the  degeneracy  of  the  bending  modes  (removing 

the  cylindrical  symmetry) ,  V  can  be  approximately 

xy 

separated  as  the  sum  of  harmonic  functions  with  different 

force  constants  (V  ~  1/2*K  (R)x^  +  1/2*K  (R)y2),  and 

xy  x  y 


. 
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the  development  which  follows  is  easily  modified] . 

An  eigenstate  of  (2.5)  must  now  be  specified  by 
three  quantum  numbers:  £,  n,  m.  An  eigenfunction  is: 

♦  £ralllr;R)  =  (z;R)  6nm(p;R)elm‘,,  •  (2Tt)"1/2  ,  (2.17) 

with  eigenvalue 

U^n(R)  =  a£(R)  +  bn(R)  ,  (2.18) 

where 

bn(R)  =  (n+l)h[kp(R)/m]1/2  .  (n  =  0,  1,  2,...)  (2.19) 

The  eigenvalues  {b^}  are  (n+l)-fold  degenerate, 
corresponding  to  the  possible  values  of  m: 

m  =  0,  ±2,  ±4,  . . .  ±n  (n  even) 

(2.20) 

m  =  ±1,  ±3,  ±5,  . . .  ±n  (n  odd) 

(Note  that  in  (2.19)  m  is  the  reduced  mass,  equation 
(2.3).  No  confusion  can  arise  from  this  dual  use  of  m; 
it  is  always  clear  from  the  context) . 

The  ($nm}  satisfy: 

{  (-h2/2m)  [32/3p2  +  (l/p)3/3p  -  m2/p2]  +  1/2  •  Kp  (R)  p2  }  Bnm 

=  VR)Bnm  '  <2-21> 

6  2 

and  are  related  to  the  Laguerre  polynomials. 


’ 


■ 
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The  "stretching  mode"  eigenfunctions  { ot ^ }  obey  the 
equation 

{ (-h2/2m) 32/3z2  +  V  z ( z ; R ) } a  ^  =  a^(R)a£  ,  (2.22) 

where  V  (z;R)  is  so  far  unspecified. 

2.  Couplings. 

Consider  the  effect  of  (2.17)  on  the  coupled  equations 
(2.12).  The  expansion  (2.4)  becomes: 


V?'R)  =  ,  *  '  (2'23) 

£,n,m  J 

and  an  element  of  the  matrix  P  (equation  (2.9))  is: 

i  i  00 

-ih<£ 'n 'm'  |  3/3R | £nm>  =  -ih6m ,m( 6r , n  /  a^,3/3Ra^dz 

—  00 

00 

+  V*  f0  B„1n3/3RBmpdp}  (2-24) 

Since,  however,  the  two  integrals  of  (2.24)  vanish 
for  £'  =  £,  and  n'  =  n,  respectively,  the  only  non-zero 
matrix  elements  are: 


—  i"h<  £  '  nm  |  3/ 3 R 


00 

£nm>  =  -ih  /  a^,3/3Ra^dz 

—  00 


(2.25a) 


and 

oo 

-ih<£n’m  |  3/3R  |  £nm>=  -ih  /  ^n*m9/aR!3nmp  dp  (2.25b) 

That  is:  vibrational  Born-Oppenheimer  states  of 
different  "stretching"  quantum  number  £  are  coupled  (2.25a), 
and  states  of  different  "bending"  quantum  number  n  are 


' 


. 
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coupled  (2.25b) ,  but  states  of  different  £  and  n  are  not 
coupled . 

This  is  not,  however,  the  same  as  saying  that  the 
bending  and  stretching  modes  do  not  interact. 

An  interaction  arises  because  an  element  (R) ,  of 
the  diagonal  matrix  U  of  equation  (2.12) ,  is  now  a  sum 
of  independent  contributions  from  both  proton  stretch 

and  bend  eigenvalues  (through  equation  (2.18)),  ap (R)  +  b  (R) . 

36  n 

In  the  Stepanov  approximation  3  Franck-Condon  pro¬ 
gressions  +  "jv^"  and  +  "jv  "  may  both  appear,  and 
the  spacing  in  each  progression  will  be  dependent  on  both 
£  and  n.  A  "stretching"  transition  from  a  curve  U^n  to 
U^,n  is  therefore  dependent  on  the  form  of  bn(R),  and 
similarly  "bending"  transitions  from  U^n  to  ,  depend 
on  (R) . 

It  has  been  observed  (a  review  is  given  by  Hadzi 
19 

and  Bratos  )  that  in  a  series  of  hydrogen  bonded  systems 
with  varying  bond  lengths,  the  frequencies  of  the  deform¬ 
ation  modes  increase  with  increasing  bond  strength 

(decreasing  R) ,  but  "that  the  increase  of  the  deformation 

19 

vibration  wavenumber  is  certainly  much  less"  than  the 

decrease  in  wavenumber  for  the  stretching  vibration. 

That  is,  Kp (R)  changes  much  more  slowly  with  R  than  does 

V  (z ; R) ;  hence  {b  (R) }  are  more  slowly  varying  than 
z  n 

(a^(R)},  and  the  couplings  (2.25b)  are  likely  to  be 
weaker  than  those  given  by  (2.25a). 


' 

' 
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An  accurate  solution  for  the  vibrational  dynamics  of 
an  isolated  hydrogen  bond  is  therefore  expected  to  be 
obtained  with  the  following  procedure: 

(i)  assume  the  separation  of  bending  and  stretching 
potentials  as  in  (2.16); 

(ii)  calculate  bending  eigenstates  in  the  Stepanov 

approximation  (retaining  R-dependence  of  K  , 

P 

and  hence  b  (R) ) ; 

n 

(iii)  solve  the  (z,R)  dynamics  with  inclusion  of 
non-adiabatic  coupling  terms,  in  order  to 
estimate  the  validity  of  the  Stepanov  approxima¬ 
tion; 

(iv)  for  isolated  cases  of  bending/stretching  near- 

1 5 

degeneracy,  allow  for  strong  coupling  by 
inclusion  of  any  non-separable  part  of  V(r,R). 

This  is  essentially  the  treatment  which  has  been 
applied  to  the  HF--  system,  except  that  dependence  of  the 
( z , R)  motion  on  the  bending  quantum  number  n  is  ignored. 
Inclusion  of  this  weak  dependence  merely  redefines  the 
potential  curves  for  heavy  particle  motion,  and  does  not 
influence  the  estimation  of  the  validity  of  the  Stepanov 
approximation  (with  which  this  calculation  is  mainly 
concerned) . 

3.  Equations  for  Stretching  Motion. 

With  the  Stepanov  approximation  for  the  bending 
motion  the  total  wavef unctions  have  the  form: 


■ 


' 
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(2.26) 


and  (2.1)  may  be  reduced  to  a  two-dimensional  equation  for 


{- (h2/2y) 92/9R2  -  (h2/2m) 92/9z2  +  V  (z,R) 

z 

+  b  (R)  -E  .  H  =  0 

n  J 


(2.27) 


The  effect  of  bn(R)  is  to  redefine  the  heavy  particle 
potential  curves.  The  index  n  is  suppressed  (assuming 
bR (R)  ~  constant  in  the  region  of  equilibrium  configura¬ 
tion)  in  the  following  development  of  the  equations  which 
have  been  solved  for  the  HF2~  system. 

It  is  then  of  interest  to  solve  the  vibrational 
dynamics  of  the  HF2-  system,  considering  motion  only  in 
the  ( z , R)  coordinates.  The  two-dimensional  Schroedinger 
equation  is: 


{- (h2/2y)  92/9R2  -  (h2/2m)  92/9z2  +  V(z,R) -E^  (z  ,R)  =  0 


(2.28) 


The  adiabatic  representation  then  uses 


(2.29) 


where  {0^}  satisfy 


{- (h2/2m) 92/9z2  +  V (z ;R) (z ;R)  =  U£(R)a£ 


(2.30) 


and  in  this  approximation  U£ (R)  =  a£ (R) . 


. 


' 


The  coupled  equations  which  result  from  the  substitu 
tion  of  (2.29)  into  (2.28)  are  again  those  shown  in 
equation  (2.12) ,  but  the  elements  of  the  matrix  P  are  now 
just  given  by  (2.25a): 

+  °o 

P£'£  “  “i*1  f  »  (z  ; R)  *  (3/3R)  z  *  (z ;  R)  *dz  .  (2 

The  procedure  for  the  solution  of  the  problem  now 
appears  to  be  straightforward: 

(i)  define  a  potential  V(z;R); 

(ii)  solve  equation  (2.30)  for  protonic  eigenvalues 
{ (R) }  and  eigenfunctions  (a^(z;R)}; 

(iii)  calculate  coupling  matrix  elements  {P^,^} 
using  (2.31) ; 

(iv)  solve  the  coupled  equations  (2.12)  for  vibra¬ 
tional  eigenvalues  {E^},  with  enough  terms  in 
expansion  (2.29)  to  ensure  convergence. 

Unfortunately,  any  attempt  to  carry  the  calculation 
"beyond  Stepanov "  encounters  an  immediate  difficulty: 
the  matrix  elements  P^,^(R)  do  not  vanish  as  R  -*■  00 .  This 
behaviour  does  not  result  from  any  approximation  in  the 
potential  function  or  the  computation;  it  is  inherent  in 
the  usual  adiabatic  representation,  and  it  is  the 
representation  which  must  be  modified. 


. 


' 
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D.  Modified  Adiabatic  Representation. 

1.  Introduction. 

In  the  adiabatic  representation  described  in  sections 
B  and  C,  the  total  wavef unction  i^(z,R)  is  expanded  in 
terms  of  vibrational  Born-Oppenheimer  (BO)  functions 
defined  on  surfaces  of  constant  R.  On  each  surface  the 
BO  functions  are  complete  in  z  space,  but  the  resultant 
coupled  equations  (in  R  space)  have  undesirable  properties. 
The  most  significant  defect,  from  the  point  of  view  of 
this  thesis,  is  that  the  matrix  P  does  not  vanish 
asymptotically . 

The  expansion  (2.29)  is  formally  correct,  and  the 
coupled  equations  (2.12)  are  mathematically  exact. 

However,  though  formally  true,  this  fact  is  not  of 
practical  use.  Problems  arise  because  individual  terms 
in  (2.29)  do  not  correspond  in  a  one-to-one  way  with 
exact  eigenstates  of  the  dissociated  systems  A-H  +  B, 

and  A  +  H-B.  This  point  has  been  discussed  in  detail  by 

mt,  ,  ^  ,  16,64 

Thorson  and  Delos. 

The  magnitudes  of  the  coupling  matrix  elements 
{P£l£}  depend,  through  equation  (2.31),  upon  the  rate  of 
change  of  the  states  {a£}  with  R  (holding  z  fixed) .  Now 
the  eigenfunctions  {a£},  being  generated  by  the  Hamiltonian 
for  the  reduced  mass  m,  must  asymptotically  (R  -*  ») 
approach  vibrational  eigenfunctions  of  the  diatomics  AH 


♦ 


. 

■ 
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and  HB  (apart  from  small  reduced  mass  errors,  of  order 
(m/y) ,  in  both  energies  and  wavef unctions ) .  For  the 
symmetric  case  a-H-A,  the  asymptotically  degenerate  (g,u) 
states  are  merely  linear  combinations  of  the  degenerate 
pair  of  diatomic  states.  In  either  case,  the  functions 
are  constant  in  form ,  since  the  potential  is  assumed  to 
take  the  dissociated  form: 

R  ->  oo 

V  ( z  ,  R)  - *  VA(ZA}  +  VB(ZB}  +  VAB(R)  (2.32) 

(z^  and  zB  are  defined  in  figure  2.2). 

Asymptotically,  a  state  merely  translates  with 
respect  to  CMN  (defined  in  figure  2.1),  from  which  z  is 
measured,  and  the  partial  derivative  of  this  function  with 
respect  to  R,  keeping  z  fixed,  is  obviously  not  zero,  but 
a  constant.  Simple  displacement  of  states  without 
distortion  cannot  give  physically  real  coupling.  There 
is,  therefore,  no  physical  basis  for  coupling  between 
asymptotic  states  in  any  reasonable  representation  of 
the  problem. 

In  the  adiabatic  representation  the  elements  (P^,^(R)} 
approach  constant  values,  the  magnitudes  of  which  (relative 
to  their  maxima)  are  not  at  all  negligible  (see  figures 
4.3  to  4.8  for  examples  of  {P^,^(R)}  calculated  using  a 
realistic  potential  for  HF2~) . 

This  phenomenon  (and  other  defects,  e.g.,  non-vanishing 
velocity  dependent  couplings  between  g  and  u  states  of 


■ 


. 

. 


Figure  2 . 2 

Jacobi  Coordinates  for  a  Linear  Bond 
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HD+^6'^  6^)  is  well  known  in  atomic  collision  theory  (a 
detailed  discussion  is  given  by  Thorson  and  Delos^6/^^)f 
and  is  the  result  of  translation  without  deformation  of 
the  wavef unctions  (this  translational  origin  is  demonstrated 
explicitly  for  the  HF system  in  Chapter  Four,  section 
C3).  The  phenomenon  does  not  appear  to  have  been  dis¬ 
cussed  in  the  literature  in  the  context  of  the  hydrogen 
bond . 

1 6 

Thorson  and  Delos  have  shown  that  a  better  representa¬ 
tion  may  be  constructed  if  the  basis  functions  are  defined 
on  surfaces  of  a  constant  coordinate  £,  related  to  z  and 

R  by  a  non-linear  transformation  (originally  suggested  by 
66-68 

Mittleman  ).  £  is  chosen  to  be  essentially  a  "heavy 

particle  coordinate",  which  is  uniquely  specified  by  the 
proton  position.  In  the  dissociated  situations  AH  or  HB, 
it  reduces  to  the  appropriate  heavy  particle  Jacobi 
coordinate  (to  be  described  below) . 

1  fi 

In  the  Thorson-Delos  (TD)  formulation,  a  wave- 
function  ipj  ( z,£)  now  satisfies  the  transformed  form  of 
the  Schroedinger  equation  (2.28),  and  may  be  expanded  in 
a  series  of  the  same  form  as  (2.29)  ,  with  £  replacing  R, 
since  the  set  of  BO  states  is  complete  for  any  parameter 
£  =  R.  The  set  of  coefficients  (f^U)}  now  satisfy  a 
different  set  of  coupled  differential  equations  ( not 
(2.12)),  and  it  is  the  purpose  of  this  section  to  exhibit 
these  equations. 


' 

■ 
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In  order  to  understand  the  motivation  for  the  TD 
choice  of  transformation,  it  is  necessary  to  examine  the 
three  sets  of  Jacobi  coordinates  and  their  relationships. 

2.  Coordinates. 

Figure  2.2  shows  the  three  possible  sets  of  Jacobi 

coordinates  for  a  linear  hydrogen  bond,  and  their 

associated  reduced  masses. 

It  will  also  be  convenient  to  use  the  position  (z  ) 

y 

of  H  from  the  geometric  centre  (GC)  of  the  heavy  nuclei 

A , B  (shown  in  figure  2.2a).  z  is  related  to  z  and  R  by: 

y 

Zg  =  z  -  1/2-AR  ,  (2.33) 

where 

X  =  (Ma-Mb)/(Ma+Mb)  ;  0  <  X  <  1  (2.34) 

A  is  called  the  "mass  asymmetry  parameter".  Without 
loss  of  generality  it  has  been  assumed:  MA  >_  Mg. 

Using  any  one  of  these  sets  of  coordinates,  the 
kinetic  energy  operator  contains  no  cross  terms: 

T  =  -  (R2/2y)  32/3R2  -  (Ti2/2m)  32/3z2  ;  (2.35a) 

or  T  =  - (h2/2pj) 32/3Rj2  -  (R2/ 2mJ)32/3zJ2  ,  (2.35b) 


where,  in  (2.35b),  J  may  be  A  or  B. 
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Now  the  defects  of  the  unmodified  adiabatic 
representation  originate  in  the  fact  that  asymptotically 
the  AH  or  HB  subsystems  are  merely  translating  with 
respect  to  the  other  atom  (B  or  A  respectively) .  This 
"free  particle  motion"  is  appropriately  described  using 
either  or  Rfi,  depending  on  the  position  of  H,  which 
is  specified  by  the  value  of  z.  That  is,  the  transla¬ 
tional  motion  associated  with  a  dissociated  diatomic  JH 

is  properly  described  in  the  Rj  system,  not  the  R  system 

1 6 

(this  is  the  concept  expressed  by  TD,  p.  152)  . 

A  transformation  is  therefore  required  which  gives 
(z,R  )  asymptotically,  and  essentially  (z,R)  in  the 
"molecular"  region. 

Using  coordinates  (z,Rj) ,  the  kinetic  energy  operator 
is,  however,  no  longer  diagonal: 


h2  fj-h2MH  .  3  .  3  ^ 

2\i  3R2  j  m(Mj+MH)  3Rj  3z  2m 


92  (2.36) 


/ 


where  fT  is  a  device 
J 

defined  by: 


introduced  for  convenience. 


64 


It  is 


f 


J 


/ 


J 


A 


/ 


J 


B. 


(2.37) 


(2.36)  is  derived  from  the  relationships: 


Mj  +  1/2 • (1-f jX)Mh 
Mj  +  Mh 


.  R  + 


f  JMH 

mj+mh  * 


(2.38) 


. 
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The  transformations  to  z T  are  given  by: 


Zj  =  z  -  1/2 •  ( A  +  f  ) R 


(2.39) 


It  is  well  known^'^0  that  by  using  mass  scaled 
coordinates,  the  motion  of  two  masses  in  three-dimensional 
space  may  be  expressed  as  that  of  a  single  particle  of 
unit  mass  on  a  six-dimensional  surface.  For  the  linear 
coordinates  used  here,  the  relationships  (2.38),  (2.39) 

may  be  expressed  as  rotations  (orthogonal  transformations) 
in  a  two-dimensional  space.  Thus  with  mass  scaled 
coordinates : 


~  _  ,  ,1/2 

z  =  (m)  •  z 


;  R  E  (y)1/2  R 


1/2  . 

‘j  ■  x“\r 


1/2 

z  -r  =  (mT)  J'/  “  *  zT  ;  Rj  =  (Pj)  .  Rj  / 


(2.40a) 

(2.40b) 


equations  (2.38)  to  (2.40)  give: 


■-V 


cos  a 


sm  a 


-sin  a 


cos  a 


J  J 


R 


(2.41) 


where 


tan  Oj  =  1/2 *(X  +  fj) (m/p) 


1/2 


=  fj(MH/MT)  (MjVMj) 


1/2 


;  (J  /  J) 


(2.42a) 


cos  cjj  =  (mj/m) 


1/2 


=  (vt/u)1/2[  (Mt  +  l/2-(l-fjA)MH)/(MT+MH)  ] 


J  H‘ 


(2.42b) 


X 


■ 
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sin  Oj  =  1/2-U  +  fj)  {irij/p ) 1/2 

=  fJ(uJ/m)1/2[MH/(MJ  +  MH)]  (2.42c) 

(2.41)  represents  two  possible  rotations:  the 
"molecular"  coordinates  (z,R)  may  be  transformed  into 
either  the  'A  configuration'  set  (z  ,R  ) ,  or  the  'B  con- 
figuration'  set  (z  ,R  ) .  In  particular,  the  second  line 

JD  ID 

of  (2.41)  suggests  the  definition  of  an  angle  ct(z;R)  with 
the  property  that  it  becomes  either  a  or  a  (depending 

Jri.  Lj 

on  the  value  of  z)  at  large  R,  for  then  a  "heavy  particle 

coordinate"  £(z;R)  may  be  constructed  which  asymptotically 

becomes  the  desired  coordinate  R  . 

J 

3.  The  Non-Linear  Transformation. 

The  preceding  discussion  is  intended  to  provide  an 

1 6 

insight  into  the  choice  by  TD  of  the  following  non¬ 
linear  transformation: 

5  =  R  cos  g(z;R)  +  z  sin  a(z;R)  (2.43) 

The  required  boundary  conditions  on  a  are  met  by  use 

of  a  switching  function  f(z;R)  (introduced  by  Schneiderman 

71  . 

and  Russek  ) ,  which  is  defined  only  through  its  asymptotic 

properties : 

Lim  [ f ( z ; R) ]  =  -1  =  f^  ;  (2.44a) 

(R  ->  zA  finite) 

and 


* 

. 
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Lim  [ f ( z ; R) ]  =  +1  =  f  .  (2.44b) 

(R  -*  °°,  zB  finite) 

a(z;R)  is  now  defined,  by  analogy  with  (2.42a): 

tan  a  ( z  ;  R)  =  l/2*[f(z;R)  +  X]  (m /\i)1/2  (2.45) 

E,  is  uniquely  specified,  at  all  R,  by  the  local  proton 
position.  In  particular,  at  large  R,  when  the  proton  is 
travelling  with  nucleus  J,  E,  becomes  the  appropriate 
coordinate  Rj.  For  the  symmetric  system  (A  =  0)  HF2-/ 
the  switching  function  shall  be  chosen  (Chapter  4,  section 
B5)  to  vanish  in  the  small  R  region.  E,  then  reduces  to 
R.  If  the  unsealed  coordinate  E,  is  defined  by: 

5  S  l/v.1/2  ,  (2.46) 

it  can  be  seen  that  E,  becomes  R  in  the  case  f  =  0,  A  =  0. 

Using  (2.43)  and  (2.45),  an  expansion  for  the 
unsealed  coordinate,  in  powers  of  (m/y) ,  gives: 

E,  =  R  +  (m/y)S(z;R)  +  0(m/y)2  ,  (2.47) 


where 

S  (z;R)  =  l/2-[  f  (z ;  R)  +A]  z  -  [  f  (z  ;R)  +A  ]  2  •  R/8  (2.48) 

S(z;R)  is  only  constrained  by  its  asymptotic  behaviour, 
and  since  f2  =  1  in  that  region,  (2.48)  may  be  written  in 
a  more  convenient  form: 


S  ( z  ;  R)  =  l/2-[ f  (z  ;  R)  + 


(z 


g 


is 


defined  in  equation 


X]  z  -  ( 1  —  X 2 ) R/8 
y 

(2.33) )  . 


(2.49) 
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Two  problems  remain.  Firstly,  the  Hamiltonian  must 
be  transformed  to  the  variables  (z,£),  and  then  new 
coupled  equations  must  be  derived  for  the  ^-dependent 
coefficients  of  a  BO  basis  (defined  on  surfaces  of  fixed 

£)  expansion  for  ^(z,£). 

1 6 

Thorson  and  Delos  obtained  the  transformed 
Hamiltonian,  retaining  terms  to  order  (m/y) ,  and  have 
shown  that  their  form  gives  correct  asymptotic  kinetic 
energy  operators  (equation  (2.36)),  within  errors  of 
order  (m/y)  . 


4.  Modified  Coupled  Equations. 

The  coupled  equations,  which  result  from  the  applica¬ 
tion  of  the  transformed  Hamiltonian  to  a  wavefunction 
expanded  in  the  BO  basis: 


(z;£) 


(£) (z; £) 


are  shown  by  TD  to  be : 


(2.50) 


{  (2y)  1  [-ifi  (d/d£)  +  P(?)  +  A(£)]2  +  U(5)}f.  (5) 

+  {A  (5)  -  (2y) -1A2  (5)  If  j  (€)  +  (2p)-1(m/y)D(S)  (-itid/dS)  2f  j  (5) 

=  E.f. (5)  (2.51) 


The  functions  -[ ct ^ }  of  equation  (2.50),  and  the 
elements  { }  of  the  matrix  U  of  equation  (2.51),  satisfy 
equation  (2.30)  with  the  replacement  R  =  £  being  understood. 


* 


Now,  in  the  final  coupled  differential  equations 
(2.51),  the  symbol  is  merely  a  dummy  variable,  and 

shall  be  replaced  by  the  more  familiar  'R',  in  order  to 
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clarify  a  comparison  with  the  unmodified  equations  (2.12). 
Evidently  there  are  additional  terms  arising  from  the 
matrices  A,  A_,  and  D.  With  the  'R'  notation,  these  have 
the  following  definitions: 

A£n(R)  =  (im/h)[U£(R)  -  Un(R)]  <  £|s|n>  ;  (2.52) 

D£n(R)  =  1/4 • <  & | f 2-l | n>  ;  (2.53) 

An  (R)  =  <Z I  V [z ; R- (m/y) S] -V (z ;R)  I n>  ;  (2.54) 

Zn  1 

1  fi  fi  4 

TD  '  have  shown  that  the  dominant  modification  to 
the  coupled  equations  is  in  the  replacement  of  the  matrix 
P  by  (P  +  A)  in  the  'momentum  operator'.  The  smaller 
terms  are: 

(i)  { A (R) - (2y ) _1A2 (R) }f j (R) ,  giving  reduced  mass 

corrections  to  the  proton  binding  energies 
(~  (m/ii)  UR)  ; 

(ii)  (2y)_1 (m/y)D(R) (-ihd/dR) 2f j (R) ,  the  diagonal 

elements  of  which  give  a  short  range  correction 
(~(m/p)E)  to  the  'heavy  particle'  kinetic  energy 
(D  vanishes  as  R  increases);  the  off-diagonal 
elements  give  couplings  which  are  an  order  of 
magnitude  smaller  than  the  (P+A)  terms. 


■ 
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Neglecting  the  smaller  terms,  the  coupled  equations 
to  be  solved  are: 

{  (2y ) _1 [  (-ihd/dR)  +  P(R)  +  A  (R)  ]  2  +  U  (R)  }f  ^  (R)  =  E^_f  (R)  . 

(2.55) 

To  repeat,  terms  neglected  in  arriving  at  (2.55)  are: 

(a)  of  order  (m/y)  ,  and  higher,,  in  the  transformed 
Hamiltonian; ^ ^ 

(b)  much  less  significant  in  the  resulting  coupled 
equations  than  the  non-adiabatic  coupling  (P+A) . 

Higher  order  terms  could,  in  principle,  be  retained, 
but  the  solution  of  (2.55)  for  a  specific  hydrogen  bonded 
system  is  expected  to  give  an  accurate  indication  of  the 
effects  of  non-adiabatic  coupling  on  vibrational  line 
positions  and  intensities.  If  these  effects  were  large, 
a  higher  order  calculation  may  be  justified,  but  for  the 
case  of  the  strongly  hydrogen  bonded  system  HF2-  (Chapters 
three  to  six)  they  are  not. 

The  purpose  of  the  TD  transformation  is  to  produce 
equations  in  which  the  principal  couplings  vanish 
asymptotically.  It  will  now  be  shown  that  the  term  (P+A) 
satisfies  this  requirement. 

5.  Asymptotic  Cancellation  of  P  and  A. 

Asymptotically,  the  vibrational  BO  states  (a£(z;R)} 
are  functions  only  of  zJf  and  not  explicitly  of  R.  Thus: 


. 


B  I 

. 
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a 


R  °o 


£ 


a£ 


(2.56) 


where  z  (z,R)  is  given  by  (2.39). 

Now  consider  the  general  transformation  (z,R)  ->  (zT/R) 

U 


(-) 


g (Zj (z ,R) ,R)  = 


P)  *  (Ji)  •  p) 

\9R  /  \9zt  /  „  \9R  / 

z  _  JR 


(2.57) 


and 


g (zJfR) 


(2.58) 


where  g  is  an  arbitrary  (well  behaved)  function. 


Thus , 


a 

n 


and 


R  -»•  oo 


a 

n 


(2.59) 


(2.60) 


The  P  matrix,  defined  by  (2.31),  becomes: 


Pn  =  -ih< £  I  ( 9/9R)  n> 
£n  1  z 


i .  e .  , 


P 


£n 


R 


00 


+  1/2- i?i  (f  j+A)  <£  I  (9/9zJ)R|n> 


(2.61) 


, 
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Also,  the  A  matrix,  equation  (2.52),  has  elements: 

A£n  =  (Wn)<£|s|n>-  (U£-Un) 

=  (im/h)<£| [h,S] |n>  ,  (2.62) 


where  h(z;R)  is  the  protonic  Hamiltonian  of  (2.30). 

Thus,  [h , S]  =  -ti2/2m{ 2  *  (  —  ^  .( — \  +  ( — } 

\9z/r  \9z/r  \9z2/r 


(2.63) 


and  from  (2.48) : 

S  ~  l/2-(fj+A)z  -  l/8-(f j+A)  2r  ,  (2.64) 

giving  32S/3z2  =  0  ,  and  3S/3z  =  l/2-(fT+A) 

J 

(2.62)  then  yields: 

A(,n  R^J°  -1/2  ifi  ( f  j+A )  <  Jl  |  (3/3Zj)R|n>  .  (2.65) 

The  required  result  is  shown  by  (2.61)  and  (2.65): 


£n 


+  A 


R  -> 


£n 


0 


(2.66) 


E.  Summary. 

To  a  good  approximation,  vibrational  eigenvalues 
{ E j }  and  eigenvectors  {fj}  satisfy  (2.55),  which  may  be 
written  in  real  form  as: 


{- (H2/2y)  [ (d/dR) I+C (R) ] 2  +  U (R) -E ^ I }f ^  (R)  =  O 


(2.67) 


where  C  is  the  real  matrix  defined  by: 


cln(R)  =i/MP,n  +  A,n) 


(2.68) 


/ 

' 
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For  a  symmetric  system  (A  =  0)  ,  A  (defined  in 

)v  n 

general  by  (2.52))  takes  the  simple  form: 

A£n(R)  =  (im/ft) (U£  -  Un)<£|0.5f (z;R) -z|n>  (2.69) 

P^n  is  defined  by  (2.31),  and  the  eigenstates  {U^; 
a^}  satisfy  (2.30)  . 

The  switching  function  f(z;R)  has  not  yet  been  defined. 
It  is  not  uniquely  specified  by  the  theory:  providing  it 
meets  the  boundary  conditions  (2.44)  it  may  be  chosen 
arbitrarily,  and  the  coupling  matrix  C  will  have  correct 
asymptotic  behaviour.  The  choice  of  a  switching  function 
for  the  HF2~  system  is  discussed  in  Chapter  4,  section  B5 . 

In  this  case  the  function  is  chosen  to  vanish  for  R  values 
at  which  the  potential  V(z,R)  has  a  single  minimum  (i.e., 
in  the  "molecular"  region,  for  which  appropriate  coor¬ 
dinates  are  z  and  R) ,  and  to  become  essentially  a  "step 
function"  (in  z  space)  as  the  barrier  to  protonic  motion 
becomes  significant. 

The  procedure  for  the  solution  of  the  vibrational 
problem  which  was  outlined  in  section  C3  (page  )  may 
now  be  adopted,  with  the  following  modifications:  in 
item  (iii)  the  integrals  {C^,^}  defined  by  equation  (2.68) 
are  required,  and  in  item  (iv)  the  coupled  equations 
(2.67)  are  to  be  solved. 

The  remaining  chapters  of  this  thesis  describe  a 
solution  of  this  system  of  equations  for  the  case  of  the 
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isolated  bifluoride  ion,  for  which  a  reasonable  potential 
function  (Chapter  Three)  has  been  constructed  from  purely 
ab  initio  data  in  the  "molecular"  region,  and  empirical 
data  in  the  dissociated  situation  (F  +  HF) .  That  the 
solution  is  convergent  is  demonstrated  by  increasing  the 
number  of  vibrational  BO  states  in  the  expansion  (2.50). 


r-  :• 


CHAPTER  THREE 


BIFLUORIDE  ION  POTENTIAL  FUNCTION 

A.  Ab  initio  and  empirical  results. 

Early  (1957-62)  investigations  of  the  electronic 

72  73 

structure  of  the  bifluoride  ion  include  valence  bond  ' 
and  SCf"^  ^  (Slater  basis)  estimates  of  the  total  energy 

O 

for  a  fixed  F-F  distance  (2.26  A),  but  no  attempt  was 

made  to  calculate  the  H-bond  energy.  McLean  and 
77 

Yoshimine  performed  SCF  (Slater  basis)  calculations 
for  21  geometries  of  HF2-/  and  found  that  the  minimum 
energy  occured  with  the  H  atom  at  the  centre  of  an  F-F 

o 

bond  of  length  2.25  A.  They  also  did  not  estimate  the 
hydrogen  bond  energy. 

24  78  79 

Later  SCF  calculations  '  '  (Gaussian  basis)  have 

shown  good  agreement  in  the  general  form  of  the  bifluoride 
ion  potential  surface:  the  minimum  energy  occurs  for  a 
linear  configuration,  with  the  proton  at  the  centre  of 

O 

the  F-F  bond,  and  interfluorine  distance  2.25-2.285  A 

°  28  80—83 

(experimental  values  for  salts  of  HF2~  are  2.26-2.29  A).  ' 

At  this  F-F  distance,  the  potential  function  for  the  proton 

is  a  fairly  flat-bottomed,  single  minimum  well  (in  good 

20 

agreement  with  the  assessment  by  Ibers  of  the  diffraction 

84 

and  spectroscopic  data;  Williams  and  Schneemeyer  have, 
however,  found  that  in  p-toluidinium  bifluoride  the  proton 

o 

is  asymmetrically  placed,  with  F-H  distances  of  1.025  A 


50 


' 
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and  1.235  A,  but  have  attributed  this  to  the  asymmetry  of 
the  crystal  environment).  As  the  F-F  distance  increases, 
a  barrier  grows  at  the  center  of  the  bond,  and  a  symmetric 
double  minimum  potential  is  formed.  Hydrogen  bond  energy 
estimates  range  from  40  to  52  kcal/mole.  Configuration 
interaction  has  not  been  included  in  any  of  these  calcula¬ 
tions  . 

At  the  time  of  writing,  there  have  been  no  entirely 
reliable  experimental  estimates  of  the  hydrogen  bond 
energy.  In  order  to  obtain  the  enthalpy  change  for  the 
gas  phase  reaction: 

HF(g)  +  F - ^HF2-(g)  (AH]_)  , 

an  estimate  of  the  lattice  energy  for  the  salt  (MHF2)  has 

8  5 

been  required  (Tuck  has  reviewed  the  situation  to  1967). 

Waddington^  calculated  U(MHF2)  (M  =  K,  Rb,  Cs)  , 

assuming  a  solid  composed  of  ionic  spheres,  and  concluded 

8  7 

AH^  =  -58  ±  5  kcal/mole.  Harrell  and  McDaniel  determined 
the  enthalpy  of  reaction  for: 

(CH3)4NF(s)  +  HF  (g)  - ►  (CH3) 4NHF2 (s)  (AH2). 

They  found  AH2  =  -37  kcal/mole,  and  assumed  the  change  in 
lattice  energy  would  be  small,  concluding  AH^  =  AH2  ±  2 
kcal/mole . 

o  p 

Dixon  et  al .  then  repeated  lattice  energy  calcula¬ 
tions  for  a  series  of  metal  salts  (M  =  Li,Na,K,Rb,Cs)  and 
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found  an  average  AH-^  =  60.2  kcal/mole.  Neckel  et  al . 

did  a  similar  calculation  (M  =  K,Rb,Cs) ,  but  used  the 

.  7  5 

wavefunction  of  dementi  and  McLean  to  calculate  the 

quadrupole  moment,  and  found  AH^  ^56  kcal/mole. 

The  "experimental"  bond  energy  estimates  for  HF2- 
are  thus  about  as  theoretical  as  those  obtained  directly 
from  ab  initio  calculation,  and  have  about  the  same  range 
of  values  (37-60  vs.  40-52  kcals/mole) . 

The  general  description  of  the  two-dimensional 
potential  surface  which  has  emerged  from  the  ab  initio 
calculations  is,  however,  intuitively  very  satisfactory: 
for  a  general  bond  A-H. . .A,  at  large  A-A  distance,  the 
proton  may  be  visualized  as  sitting  in  either  of  two 
identical  Morse-like  curves;  this  is  a  double  minimum 
potential,  the  barrier  height  being  the  dissociation 
energy  De,  for  the  diatomic  A-H.  As  the  two  A  nuclei 
approach,  the  barrier  height  must  decrease  to  zero  (al¬ 
though,  of  course,  the  minimum  energy  does  not  necessarily 
occur  when  there  is  no  barrier) ,  for  at  very  small  A-A 
distance,  a  barrier  to  protonic  motion  is  physically 

inconceivable . 

2  0 

Ibers  proposed  an  empirical  form  for  the  potential 

function  governing  the  stretching  motions  of  the  bifluoride 

ion.  In  his  polynomial  representation,  the  coupling  of 

the  two  modes  (Q-^  for  symmetric,  for  asymmetric,  normal 

2 

coordinates)  is  included  through  a  term  dQ^Q^: 
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V(Q1,Q3)  =  aQ^  +  bQj  +  cQ3  +  dQ1Q3  (3.1) 

Ibers  obtained  solutions  to  the  Schroedinger  equation 
with  this  potential,  by  using  an  expansion  of  products  of 
harmonic  oscillator  functions  (Chapter  One,  section  B2); 
the  four  parameters  (a,b,c,d)  were  varied  to  give  agree¬ 
ment  with  the  following  four  experimental  values: 

(symmetric  stretch,  600  cm  (asymmetric  stetch, 

1450  cm  ;  a  value  assigned  to  the  second  overtone,  3v^ 

(5090  cm  ■*■)  ;  and  an  estimation  of  the  contraction  A,  in 

O 

the  equilibrium  F-F  distance  upon  deuteration  (0.0024  A). 

A  computational  error  in  the  first  estimation  of  these 

'force  constants'  was  later  corrected  by  Ibers  and 
27 

Delaplane . 

7  8  7  9 

Kollman  and  Allen,  and  Noble  and  Kortzeborn  have 

attempted  to  fit  their  ab  initio  data  points,  and  those 

77 

of  McLean  and  Yoshimine,  to  the  form  of  Ibers.  The 

resultant  force  constants  are  quite  dissimilar.  Noble  and 

Kortzeborn  have  obtained  coefficients  for  a  higher  order 

3  4 

polynomial  representation  (including  and  Q^) ,  and  have 
compared  these  coefficients  with  those  obtained  using  the 
data  of  McLean  and  Yoshimine.  Again  the  agreement  is 
very  poor,  and  the  correspondence  between  coefficients 
of  different  polynomial  representations  is  even  worse. 

Janoschek  has  summarized  the  situation:  "the 
potential  parameters  for  different  analytical  representations 


I 
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of  an  energy  surface  are  only  of  little  importance... 
furthermore  the  same  analytical  representations  are  unable 
to  be  compared...  only  the  calculated  frequencies  are 
able  to  yield  a  global  judgement  for  the  accuracy  of  the 
representation. " 

Alml0f^4  has  performed  a  calculation  of  the  potential 
surface  and  the  associated  vibrational  levels.  His 
MO-LCAO-SCF  (Gaussian  basis)  calculation  is  the  most 
extensive  of  those  discussed,  involving  27  linear  geometrie 
of  the  bifluoride  ion.  Again,  a  polynomial  representation 
for  the  vibrational  potential  surface  was  used,  and  the 
coefficients  were  obtained  by  least-squares  fit  of  his 
ab  initio  points.  The  minimum  of  this  surface  occurs  at 

o 

an  F-F  distance  of  2.246  A,  in  a  single  minimum  region. 

The  vibrational  problem  was  also  solved  by  expanding  the 
wavefunction  in  products  of  linear  ho  functions,  and 
convergence  to  ±5  cm  ^  was  obtained  with  160  products  (!). 
The  frequencies  reported  by  Alml0f  are  in  good  agreement 
with  experimental  data.  For  the  lower  energy  transitions, 
at  least,  the  polynomial  representation  appears  to  be  a 
reasonable  potential  surface.  No  calculation  of  relative 
intensities  was  performed. 

For  the  purpose  of  the  calculation  described  in  this 
thesis,  all  polynomial  representations  of  the  function 
V ( z , R)  have  a  serious  fundamental  defect:  they  take  no 
account  of  dissociation  in  the  coordinate  R.  Furthermore, 
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beyond  a  fairly  small  region  of  R  (in  which  they  are,  in 
various  ways,  constrained  to  be  reasonable  representations), 
they  may  exhibit  wildly  non-physical  behaviour  (polynomial 
oscillations) ,  for  which  the  vibrational  Born-Oppenheimer 
eigenvalues,  and  matrix  elements  {  (i/ft)  P^n  (R)  } ,  will  be 
completely  unreliable. 

An  analytic  form  for  the  potential  function  of  a 
general  hydrogen  bond,  involving  the  superposition  of 

Morse-like  curves,  was  proposed  by  Lippincott  and 
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Schroeder,  and  modified  by  later  workers.  '  This 

model  is  purely  empirical,  requiring  information  from 
experiment,  and/or  ab  initio  calculation,  in  order  to 
establish  the  values  of  several  parameters  for  a  particular 
bond.  It  does,  however,  have  the  essential  characteristics 
of  reasonable  form  in  the  dissociation  limit,  and  decreas¬ 
ing  barrier  height  in  the  interaction  region.  The  form 
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(but  not  the  parameters)  of  Jiang  and  Anderson,  for  this 
representation  of  the  H-bond  potential  function,  has  been 
adopted  for  the  calculation  described  herein. 

B.  An  Asymptotically  Reasonable  Form. 

24 

1.  The  Polynomial  Form  of  Alml/f. 

Before  discussing  the  final  form  for  the  potential 
used  in  this  study,  it  is  necessary  to  describe  the 
polynomial  obtained  by  Alml/f,  and  its  characteristics, 
because  the  free  parameters  in  the  Jiang  and  Anderson  form 
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have  been  fixed  using  this  ab  initio  data. 

Figure  3.1  shows  the  relationships  between  the 

coordinates  used  by  Alml/f  (Q-^C^)/  Jiang  and  Anderson 

(s,t) ,  and  this  work  (z,R). 

Alml^f's  polynomial  is  a  more  general  form  of  that 
2  0 

of  both  Ibers  (equation  (3.1)),  and  Noble  and 
79 

Kortzeborn.  It  is: 


(3.2) 


(Note  that  here  Q2  is  used  in  place  of  for  the 
asymmetric  coordinate) . 

The  paper  by  Alml/f  contains  several  errors:  in  his 
definition  of  Q2 ,  the  factor  of  1/2  is  missing  (see 
figure  3.1);  in  'figure  1'  of  Alml/f's  paper,  both  axes 
are  incorrectly  labelled  (this  is  a  contour  diagram  for 
the  potential  surface;  the  axis  labelled  'Q^(a.u.)',  should 


o 


be  'Q^  +  Re  (A)',  and  the  other,  labelled  'Q2(a.u.)', 


O 


should  be  ' 2Q2  (A)').  These  errors  have  been  resolved  by 
private  communication  with  Dr.  Alml/f. 

It  is  convenient  to  express  V(Q1,Q2)  in  a  form  which 
emphasizes  the  non-separable  part.  The  full  expansion  is: 


V(Ql'Q2)  Ql^C20  +  C30Q1) 


(3.3) 


Now  defining: 


Q  E  Qx  =  R-R 


(3.4) 


e 
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Figure  3.1 

Coordinates  for  the  Linear  Bifluoride  Ion 


2 


Centre  of  Mass  of  F^,  F2* 

The  subscript  'e'  is  used  to  denote  'equilibrium 
values',  that  is,  values  for  the  coordinates  at  the 
minimum  of  the  potential  surface. 

Thus : 

R  =  2 . 246  A  (4 . 244  a.u. ) 
e 

s  =t  =  1.123  A 
e  e 

z  =0.0 
e 

Alml^f  coordinates: 

Qx  =  R-Re  =  (s-se)  +  (t-te) 

Q2  =  (s-t)/2  =  z. 

Jiang  and  Anderson  coordinates: 
s ( z , R)  =  R/2  +  z 
t ( z , R)  =  R/2  -  z 


■R 


R/ 


Fi 


H 

i 


. 
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(3.3)  may  be  written: 

V(Q,z)  =  V^(R)  +  V^(z,R) 


(3.5) 


with: 


V^(R)  =  Q2(c20  +  C30*Q)  (3.6a) 

V^(z,R)  E  Az6  +  B  (Q)  z 4  +  C(Q)z2  (3.6b) 


and : 


B(Q)  5  CQ4  +  C14-Q 

(3.7) 

C(Q)  E  CQ2  +  Cl2Q  +  C22Q2 

Clearly,  V^(z,R)  is  symmetric  in  the  z  coordinate, 
at  any  R.  It  is  the  variation  in  the  coefficients  of 
this  sixth  degree  polynomial  which  causes  the  growth  of 
a  barrier  as  R  increases.  The  partial  derivative, 

3V^/3  R,  will  influence  the  form  of  the  matrix  elements 
{  (i/tiJP^tR)  }. 

Figures  3.2  to  3.4  show  the  functions  (R) ,  and 
V^(z,R)  over  a  range  of  values  of  R. 

In  the  region  R  =  3.8  to  5.2  a.u.,  where  Alml^f  has 
ab  initio  data  points,  the  polynomials  V^,  and  have 
reasonable  physical  behaviour.  Beyond  R  =  5.2  a.u.  the 
cubic  behaviour  of  is  evident  in  figure  3.2,  and  the 
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FIGURE  3.2 


V^(R)  x  10^  ( au. )  vs  R  (a.u.) 
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FIGURE  3.3 

yA(z;R)  x  103(a.u.)  at  R  =  4  to  5a.u. 
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FIGURE  3.4 
A  3 

Vj(z;R)  X  10  (q.u.)  at  R  =  5.6,  5.8  a.u. 
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oscillations  of  the  sixth  degree  polynomial  V^,  become 
apparent  in  figure  3.4. 

To  summarize,  the  Alml^f  polynomial  representation 
is  physically  reasonable  in  a  small  region  of  R,  about 
the  equilibrium  position  (R  =  4.244  a.u.),  but  has 
unacceptable  asymptotic  (dissociation  in  R)  character. 

That  Alml^f's  calculated  vibrational  frequencies  are 
accurate,  for  lower  energies  (within  5%  of  experimental 
values) ,  is  taken  to  indicate  that  the  surface  is  a  good 
representation  in  this  region. 
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2.  The  Jiang  and  Anderson  Form. 

90  91 

Lippincott  and  Schroeder  '  have  shown  that  the 
vibrational  potential  function  for  a  diatomic  molecule 
(AB)  is  well  represented  by  the  Morse-like  form: 

V7,D(r)  =  De  (1-exp  (-nAr2/2r)  )  (3.8) 

Ad 

where  r  is  the  internuclear  distance,  Ar  =  r-re  is  the 
displacement  from  equilibrium  (VAB(re)  E  0)  »  an<^  n 
obtained  from  the  curvature  at  re  (empirically,  using  the 
force  constant  k^B  from  the  IR  fundamental,  n  =  k^B*re/De) . 

Attempts  were  made  '  '  to  simulate  the  potential 
function  for  a  linear  H-bond,  A-H. .A,  by  superposition  of 
forms  (3.8),  together  with  terms  intended  to  include 
electrostatic  attraction  (-C2 (Re/R) m) /  and  Van  der  Waals 
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repulsion  (C^exp (-bR) ) .  Thus  the  Lippincott  and  Schroeder 


■ 
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(LS)  form  for  the  H-bond  potential  is: 

VLS  =  0(1  "  y(s)  “  y(t)  +  ~  C2(Re/R)m)  (3.9) 


where 


Y  (u)  =  exp(-n(u-ueJ  /2u)  (3.10) 

ueoo  =  equilibrium  bond  length  for  free  diatomic  A-H. 

(seoo  =  teoo  =  °-917  A  for  HF) 

n  =  k„_-u  /D 
HA  e°° 

(kHF  =  9.655  mdyn.  A  ^ ,  D  =  51,300  cm  b , 
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dissociation  energy  for  H-F  bond) 


Re  =  equilibrium  H-bond  length. 

O 

(R  =  2.246  A,  from  Alml^f,  has  been  used) 

C^,  C 2,  b  and  m  must  also  be  fixed  empirically. 

In  attempting  to  fit  the  ab  initio  data  of  Kollman 
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and  Allen,  and  Alml^f,  Jiang  and  Anderson  have 

found  the  form  (3.9)  to  be  insufficiently  flexible,  and 

have  introduced  an  additional  repulsive  term  g-y(s)-y(t): 


vja(2'r)  =  D(l-y (s)-y  (t)+gy (s)y (t)+C1e  bR  -  C2(Re/R)m) 

(3.11) 

There  are  now  five  undetermined  parameters:  C-^,  C 2,  g, 
b,  and  m. 

C^,  and  C2  may  be  fixed  from  the  following  equilibrium 
conditions : 


V  (0,R  )  =  Ve  =  hydrogen  bond  energy 
u  A  0 


(3.12a) 


(-51  kcal/mole,  from  Alml^f,  was  used) 


■ 


and , 


0 


(3.12b) 


8VT7V 
- JA 

8s 


se ' 


3V 

3t 


JA 


s  ,  t 
e  e 


giving : 

C2  =  (Ve/D-l-y (se)  (g*y  (se)-2)-y’  (se)  (gy (sQ) -l)/b) / (m/bRe-l) 

(3.13) 

hR 

C-l  =  e  e(mC2/Re+y'  (se)  (g*y(se)-l)  )/b  (3.14) 

Thus,  given  the  parameters  m,  g,  and  b,  the  equilibrium 
conditions  (3.12)  have  been  used  to  determine  the  coef¬ 
ficients  C^,  and  C2 • 

Jiang  and  Anderson  have  assigned  values  (m  =  6 , 
g  =  0.64,  b  =  3.2009  A  ^)  by  fitting  the  form  (3.11)  to 
ab  initio  data,^'^  and  "observed  infrared  frequency". 

In  this  study >  the  surface  of  Almldf  alone  has  been 
used  to  define  the  parameters. 

First,  the  parameter  g  is  fixed  using  V^(z,R)  (equation 
(3.6b),  and  figure  3.3).  For  comparison,  V^(z,R)  is  defined 
by: 


V1(z,R)  =  Vja(z,R)  -  Vja(0,R)  (3.16) 

This  merely  sets  V1(0,R)  =  0,  for  all  R  (as  V^(z,R)),  and 
using  (3.11)  this  becomes: 


„ 
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V1(z,R)  =  D (g  (y (s) y (t) -y2 (C) ) +2y (O -y (s) -y (t) )  (3.17) 

where , 

£,  =  s  (0 ,  R)  =  t  (0 ,  R)  =  R/2  (3.18) 

Hence,  V^(z,R)  is  only  dependent  on  g  (not  on  m  or 
b) .  With  a  value  g  =  0.55,  equation  (3.17)  has  been  found 
to  give  curves  which  are  in  good  agreement  with  those  of 
figure  3.3,  in  the  region  R  =  3. 8-4. 8  a.u. 

Now  V 2 (R)  is  also  defined  to  be  consistent  with  the 
Alml^f  function  (equation  (3.6a),  and  figure  3.2)  V^R): 

VR)  5  VJA(0'R)  -  VJA(0'Re>  (3-19) 

so  that  V0(  R  )  =  0 . 

2V  e 

The  complete  potential  surface  is  now: 

V (z , R)  =  V2(R)  +  Vx(z,R)  (3.20a) 

i  .  e .  , 

V (z ,R)  =  Vja(z,R)  -  VJA(0,Re)  (3.20b) 

(3.20b)  simply  defines  the  potential  to  be  zero  at 
the  equilibrium  position  for  the  bifluoride  ion  (as 
Alml^f ) ,  rather  than  at  an  equilibrium  configuration  for 
the  dissociated  fragments,  F-  +  HF ,  or  FH  +  F-  (as  Jiang 


and  Anderson) . 


Having  fixed  g,  the  function  V^(z,R)  is  fixed,  and 

the  vibrational  Born-Oppenheimer  eigenvalues  (e^(R)}  may 

be  calculated  by  the  method  of  Chapter  four.  The  curves 

(U^(R)}  may  then  be  constructed  (U^ (R)  =  (R)  +  V2 (R) ) 

for  a  given  pair  of  the  parameters  m  and  b.  The  values 

of  these  parameters  were  varied  (only  integer  values  of 

m  were  considered)  until  the  curves  (U^(R)}  agreed  closely 

with  those  calculated  using  the  Alml^f  surface  in  the 

physically  reasonable  region,  R  =  3. 8-5. 2  a.u.  In  this 

°  _1 

way,  the  values  m  =  4,  b  =  2.35  A  ,  have  been  assigned. 

3.  Summary. 

A  potential  function  V(z,R)  has  been  defined  (equation 

(3.10),  (3.11)  and  (3.20b))  using  the  representation  of 
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Jiang  and  Anderson,  with  parameters: 

m  =  4 

°-l 

b  =  2.35  A 
g  =  0.55. 

This  representation  gives  a  potential  surface  which 

2  4 

matches  that  of  the  polynomial  form  of  Alml^f  very  well 
in  the  region  3. 8-5. 2  a.u.,  and  in  addition  gives  good 
behaviour  as  R  00 . 

It  should  be  noted,  however,  that  at  small  R3  the 

function  is  dominated  by  the  attractive  term  -C~  (R /R)m 

JA  ^ 


Thus , 


. 


1 1 


67. 


Lim  V 2  (R)  =  -  00  , 

R+0 

and  the  H-bond  collapses.  The  Jiang  and  Anderson  representa¬ 
tion  breaks  down  at  small  R.  This  does  not,  however, 
present  any  difficulty  in  practice,  because  the  matrix 
elements  { (i/ft) Pmn (R) } ,  and  curves  (Un(R)}  may  be  numerically 
extrapolated  without  difficulty  for  R  <  3.8  a.u.,  in  a 
manner  which  avoids  this  false  behaviour  (Chapter  5, 
section  G) . 


CHAPTER  FOUR 


VIBRATIONAL  BORN-OPPENHEIMER  STATES 

A.  General  Eigenvalue  Condition. 

This  chapter  describes  a  method  for  the  numerical 
solution  of  the  eigenvalue  problem  (at  fixed  R)  posed  by 
equation  (2.30),  with  V(z;R)  =  V][(z;R)  +  V2  (R)  ,  and 
(R)  =  £^(R)  +  V2 (R) ,  as  described  in  Chapter  3. 

For  an  arbitrary  energy  e,  there  exist  solutions  to 
the  differential  equation: 

{3^/9z^  +  (2m/h^) [e-V^ (z ;R) ] }a (z ;R)  =0  .  (4.1) 

Unless  e  is  an  eigenvalue,  however,  a(z;R)  will  not 
be  everywhere  'regular',  that  is,  will  not  satisfy  the 
boundary  conditions: 

a(z;R)  Z  — >  0  (4.2) 

An  eigenfunction  a^(z;R)  is  everywhere  continuous, 
obeys  (4.2),  and  satisfies  (4.1),  for  which  the  correspond¬ 
ing  energy  is  an  eigenvalue  e^(R).  R  is  here  considered 
to  be  a  fixed  parameter  (indicated  by  the  semi-colon) , 
and  so  the  quantum  states  {a^}  °f  the  reduced  mass  m  are 
essentially  vibrational  Born-Oppenheimer  states  of  the 
proton. 

The  JWKB^^  approximation  guarantees  that,  in  classically 
forbidden  regions  of  z  far  from  the  turning  points ,  the 
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form  of  a(z;R)  for  an  arbitrary  energy  e,  is: 


(4.3a) 


where  c  is  an  arbitrary  constant,  and 


p(z;R)  5  {  (2m/fi2)  [V1(z;R)-e]  }1/2 


(4.3b) 


There  is,  therefore,  one  solution  which  is  exponentially 
growing  and  one  which  is  exponentially  decreasing  in  a 
direction  away  from  either  turning  point.  Growing  solutions 
are  clearly  forbidden  by  (4.2). 

At  any  point  in  z  space,  z^,  a  solution  to  (4.1)  is 
uniquely  defined  by  the  values: 

ou  =  a (z.;R) , 

and  a  .  ?  =  a  '  (z  .  ;  R)  .  (4.4) 

l  l 

If  z^  is  in  the  right  asymptotic  region  (RAR:  z-H-°°)  , 
and  is  any  arbitrary  value,  then  the  derivative  is 
related  to  this  by  (4.3): 


(4.5a) 


Similarly,  in  the  left  asymptotic  region  (LAR:  z+-°°)  : 


(4.5b) 


In  obtaining  (4.5)  from  (4.3),  the  JWKB  criterion 
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(4.6) 
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has  been  used. 

The  different  signs  in  (4.5a)  and  (4.5b)  ensure  the 
correct  asymptotic  exponential  decay  of  the  solution. 

Solutions  to  (4.1)  may  then  be  accurately  numerically 
propagated  from  either  of  the  points  or  zv'  towards 
the  left  and  right  turning  points  respectively.  In  fact, 
the  initial  values  assigned  to  and  cad  are  quite 
unimportant  (providing  they  are  not  both  zero)  because 
any  non-zero  pair  of  numbers  may  be  considered  to  be  a 
linear  combination  of  the  two  solution  forms.  Propagation 
towards  the  turning  points  ensures  the  decay  of  the 
unacceptable  form,  and  the  process  is  self-purifying. 

Thus,  if  solutions  a  and  a  are  generated  regular 
in  LAR  and  RAR  respectively,  they  may  be  accurately 
propagated  to  a  suitable  matching  point  z  .  zm  must  be 
chosen  in  the  classically  allowed  region,  because  a  solu¬ 
tion  which  is  regular  in  one  region  cannot  be  propagated 
into  the  other  classically  forbidden  region  without  the 
undesirable  exponential  growth  occurring. 

The  condition  that  e  be  an  eigenvalue  is  clearly 

that  the  two  solutions,  generated  from  LAR  and  RAR,  be 

identical,  for  then,  and  only  then,  is  the  function 

regular  on  —  °°  to  +°°.  Thus,  functions  and  derivatives 

must  match  at  z  . 

m 

The  eigenvalue  condition  is: 


' 


‘ 
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aaL  +  baR  =  0  ; 

L  •  pi 

and  aa  +  ba  =  0  (4.7) 

That  is,  the  constants  a  and  b  must  exist,  which 
requires : 


(4.8) 


This  is  the  well  known  property  of  the  Wronskian  of 

two  linearly  dependent  solutions  to  a  second  order  linear 

differential  equation:  the  Wronskian  vanishes  at  any 
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point . 

If  the  function  W(e)  is  defined: 

,,/v—LR1  R  L '  //in\ 

VI  (e)  -  cl  a  -  a  a  (4.9) 


the  eigenvalue  condition  becomes: 

W(e£)  =  0  (4.10) 

W(e)  is  to  be  calculated  at  the  matching  point  z^. 

These  concepts  will  be  extended,  in  Chapter  5,  for 
the  case  of  N  coupled  second  order  differential  equations. 


(B) .  Symmetric  Potential: 

1.  Constraints: 

V1(z;R)  has  been  defined  such  that: 
V1(0;R)  =  0  ; 
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and 


V1(z;R)  =  V1(-z;R)  . 


(4.11) 


This  leads  to  the  eigenfunction  conditions: 
Even  parity: 


(0 ; R)  =  0 
(0; R)  =  C 


U  =  0,2,4...) 


(4.12a) 


Odd  parity: 


a ^ (0;  R) 
a£'  (0;  R) 


d 


0 


(£  =  1,3,5...) 


(4.12b) 


c  and  d  are  arbitrary  initial  values  (eventually  fixed 
by  normalization) . 

2.  Numerical  Method: 

The  conditions  (4.12)  lead  to  a  great  reduction  of 
numerical  computation,  particularly  for  cases  in  which 
V-^(z;R)  has  a  single  minimum,  or  the  top  of  the  barrier  is 
below  the  eigenvalue  being  calculated. 

Numerical  calculation  of  an  eigenvalue  falls,  then, 
into  two  classes. 

(i)  e  _>  0 : 

In  this  case  the  matching  point  is  chosen  at  z  =  0 . 
It  is  only  necessary  to  propagate  a  solution  from  one  of 
the  asymptotic  regions,  the  other  solution  being  fixed 
by  (4.12). 

Arbitrarily  choosing  to  propagate  from  RAR,  (4.9) 


becomes : 
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D( 

W(e)  =  a  even  parity 

W(e)  =  aR  odd  parity  (4.13) 

Values  of  e  are  sought  such  that  either  the  wave- 
function  vanishes  at  z  =  0  (odd  parity  eigenvalues) ,  or 
the  derivative  vanishes  there  (even  parity  eigenvalues). 

The  fourth-order  Runge-Kutta  procedure  has  been  used 
to  integrate  (4.1)  from  RAR  to  z  =  0 .  The  Newton-Raphson 
method  was  used  to  find  the  zeros  of  W(e). 

Eigenvalues  have  been  calculated  to  eight  significant 
figures  by  this  method,  and  are  completely  stable  on 
halving  the  Runge-Kutta  integration  stepsize  and  varying 
the  starting  point  in  RAR. 

(ii)  e  <  0 : 

z  =  0  is  not  a  suitable  matching  point  if  the  distance 
of  the  energy  from  the  minimum  of  the  potential  well  is 
less  than  the  barrier  height.  This  is  always  the  case  for 
bound  states  at  sufficiently  large  R.  The  problem  is 
that,  in  order  to  reach  z  =  0,  the  solution  propagated 
from  RAR  must  penetrate  a  classically  forbidden  region  in 
a  direction  away  from  a  turning  point,  and  then  exponential 
growth  of  the  forbidden  solution  occurs.  This  is  overcome 
by  propagating  a  solution  from  z  =  0,  with  initial  values 
defined  by  (4.12),  towards  the  RAR  solution.  The  solutions 
are  then  matched  in  a  classically  allowed  region,  at 


*. 


' 


74  . 


z^  =  Zq  (the  position  of  the  potential  minimum) ,  and  zeros 
are  sought  for  W(e)  defined  by  the  full  equation  (4.9). 


3.  Eigenvalue  Initialization  : 

In  seeking  a  zero  of  W(e),  the  value  of  e  must  initially 
be  chosen  close  to  the  required  eigenvalue  e^,  and  must  be 
limited  by  high  and  low  bounds,  so  that  computational  effort 
is  minimized  and  convergence  to  an  undesired  eigenvalue  is 
avoided.  This  is  achieved  by  starting  the  calculation  of 
(e^(R) ;a^(z;R) }  in  a  region  of  R  where  the  potential  has  a 
single  minimum  and  non-zero  curvature  at  z  =  0 .  A  harmonic 
approximation  then  yields  a  good  approximation  to  £q,  and 
subsequent  eigenvalues  are  initialized  from  the  spacing  of 
the  previous  ones .  Providing  the  step  to  the  next  R  point  is 
sufficiently  small,  e£(Ri+1)  is  approximately  e^(Ri)  and  the 
process  continues. 


Alternatively,  the  procedure  may  be  initiated  in  a 
region  of  large  R,  where  the  double  minimum  gives 
effectively  harmonic  oscillator  states  of  the  non— hydrogen 
bonded  A-H  molecule. 


4.  Normalization; 

The  wavef unctions  (a^(z;R)}  are  stored  pointwise  in 

z  space,  and  (after  suitable  matching  of  magnitudes  at 

z  for  the  case  £0  <  0)  each  is  numerically  integrated 
m  ~ 

(Simpson's  rule  is  found  to  be  quite  adequate  for  six 
figure  accuracy  in  the  subsequently  calculated  P  and  A 
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matrix  elements)  to  find  the  scaling  factors  C^,  such  that: 

2  •  /°°  (a «  •  c71//2 )  2  dz  =  1  .  (4.14) 

0 

The  functions  are  thus  redefined  as  the  normalized 
set  {a^}. 


5 .  Matrix  elements :  P . . ;  A . . ;  C . . . 
- ij - ID- - id- 


(i)  P-  • 
13 
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The  Hellmann-Feynman  theorem  gives: 


■00 


(i/h)  *P£n(R)  *  [en(R)”^(R)]  =  -t  a£(z 


;R)(aR  )  an(z;R)dz* 


(4.15) 


53 

V^,  using  the  Jiang  and  Anderson  form,  has  been  defined 
(equation  (3.17))  by: 

V1(z;R)  =  D{g*  (y (s)y (t) -y (£) 2)+2y  (£) -y  (s) -y  (t) }  (4.16) 

with  £  =  R/2. 


Any  part  of  8V1/8R  which  is  only  a  function  of  R  gives 
zero  contribution  to  (4.15).  Now: 

[9V-./3R]  =  (D/2 )  •  d  (z ;  R)  +  function  of  R  (4.17) 

JL  Z 

where : 

d ( z ; R)  E  y '  (s) {g*y (t) -1}  +  y '  (t) {g-y  (s) -1}  (4.18) 

and  (4.15)  becomes 

(i/h)Pn  (R)  *  (e  -e  0 )  /D  =  f°°  a  n  (z ;  R)  d  (z ;  R)  a  (z  ;  R)  dz  (4.19) 
x,n  n  X/  n  “ 
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This  integral  is  calculated  numerically  at  each  R 
point  using  the  stored  normalized  {a£}.  It  is  evident 
from  (4.15)  that  there  is  no  coupling  between  states  of 
opposite  parity  ( 9V^/3R  is  an  even  function) . 

(ii)  A . . : 

13 

From  equation  (2.69),  the  matrix  elements  required 
to  correct  for  proton  translation,  are  given  by: 

+00 

(-ih)A^n(R)  =  m  (££-en)  *  1/2  *  • z • f (z ;R) a^dz  (4.20) 

The  form: 

f ( z ; R)  =  tanh (y (R) -Zq (R) • z)  (4.21) 

has  been  chosen.  This  has  the  property  that  when  there  is 
no  barrier  (z^  =  0)  the  switching  function  vanishes. 

y (R)  has  been  chosen  in  such  a  way  that  all  elements 
of  the  matrix  C(R)  remain  smooth  while  asymptotically 
going  to  zero. 

The  following  form  (shown  in  figure  4.1)  was  found  to 
be  suitable: 

y(R)  =  a{ 1  +  tanh (b (R-R  ) ) }  (4.22) 

-2  -1 

with:  a  =  1.7  a.u.  ;  b  =  5.0  a.u.  ;  Rc  =  4.97  a.u. 

It  is  not  claimed  that  the  parameters  a,  b,  Rc  are 
unique.  They  are  merely  convenient  in  defining  a  switching 
function  which  grows  rapidly  from  zero,  in  the  single 
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FIGURE  4-1 

GAMMA  (a.u.)~^  vs.  R  (a.u.) 
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minimum  molecular  situation,  to  ±1  in  the  high  barrier 
(dissociated)  situation.  The  form  (4.22)  provides  a 
continuous  transformation  between  these  two  limiting 
situations . 

From  (4.20)  the  integrals  to  be  numerically  evaluated 

are : 

(-ifi)  A^n/ (m  (e^-en)  )  =  f  a^zfo^dz  (4.23) 

(iii)  C .  .  : 

ID 

The  C  matrix  is  defined  to  be  real,  by: 

CJtn(R)  =  +  A*n(R))  (4’24) 

C.  Results. 

(en(R)}  are  eigenvalues  of  V^(z;R),  but  the  full 
potential,  V(z,R),  has  been  written  as  (z } R) +V2 (R) , 
and  so,  in  order  to  provide  potential  surfaces  for  the 
complete  eigenvalue  problem, (U^ (R) }  are  now  defined  to  be 
(e£(R)+V2(R)}. 

1.  hf2~. 

Figure  4.2  shows  the  curves  (U^(R),  %  =  0,5}  obtained 
for  HF2-  using  the  potential  function,  of  the  Jiang  and 
Anderson^  form,  described  in  Chapter  3.  Jiang  and 
Anderson  obtained  similar  curves  using  a  different  numerical 
method  for  the  solution  of  (4.1).  They  did  not  calculate 
matrix  elements  or  attempt  to  solve  the  coupled  equations 
for  vibrational  eigenvalues.  They  discussed  'vertical 
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FIGURE  4.2 


80. 


transitions'  from  the  minimum  of  each  Un (R)  curve.  This 

is  incorrect,  and  is  a  poor  approximation. 

Figures  4.3  to  4.5  show  matrix  elements  {P. .},  and 

13 

{Cij},  coupling  the  first  three  states  of  even  parity, 
while  figures  4.6  to  4.8  show  the  odd  parity  set. 

The  maximum  of  is  at  about  4.75  a.u.,  whereas 

the  minimum  of  the  curve  Uq (R)  is  at  4.30  a.u.  It  is  then 
immediately  clear  that  a  ground  state  wavefunction  in 
the  curve  Uq  will  be  only  weakly  coupled  to  any  function 
associated  with  curve  .  Similar  considerations  for 
and  indicate  that  a  ground  state  wavefunction  in 
will  probably  be  weakly  coupled  to  functions  associated 
with  . 

It  is  therefore  expected  that  the  fundamental  I.R. 
transition  (v^)  will  be  well  approximated  by  the  difference 
in  energy  between  a  harmonic  oscillator  (ho)  approximated 
ground  state  in  U^E^Oio))  and  a  similar  ground  state  in 

U0(E00(ho)): 

V3  -  Eoi^ho^  “  E00^h°^  (4.25) 

Furthermore,  the  Raman  active  transition  (v^  should 
be  well  approximated  by: 

V1  -  E10(ho)  ”  E00^h°^ 

where  ElQ(ho)  is  the  first  excited  state  in  the  ho 
approximation  for  Uq (R) . 


(4.26) 
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FIGURE  4-3 


i/frf^2  &  C  (au.)'1  vs  R(a.u.) 
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FIGURE  4-4 
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FIGURE  4-5 


(a.u.)"^  vs  R(a.u.) 
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FIGURE  4.6 

(i/till’j  &  C  (au.)"'  VS.  R(au) 


R(aa) 
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FIGURE  4.7 

-  (r/ti)  P  &  -C  (a.uj"1  VS.  Rfcuj.) 
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FIGURE  4.8 
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It  is  not  clear  that  higher  eigenstates  of  will 
have  negligible  coupling  to  states  of  U^,  for  the  higher 
wavef unctions  have  a  considerable  region  in  the  classic¬ 
ally  allowed  domain  which  overlaps  the  maxima  of  C.^,  C^,_, 
and  The  quantitative  effects  of  this  are  revealed 

by  the  accurate  solution  of  the  coupled  equations,  the 
subject  of  Chapters  5,  and  6. 

2.  Deuterium  Shifts: 

Curves  (U^(R)},  for  DF2~  are  readily  obtained  by 

solving  (4.1)  with  M  for  M  in  the  reduced  mass  m.  These 

u  H 

curves  are,  of  course,  more  closely  spaced  than  those  for 
HF^-/  (U^(R)},  but  are  otherwise  similar,  and  have  not 
been  shown  here. 

The  non-adiabatic  coupling  curves,  (P^j(R)},  have 
slightly  higher  maxima  (as  expected  from  (4.15);  the 
maximum  of  (i/hjP^/  for  example,  is  about  4.16  (a.u.)  ^)  , 
but  are  otherwise  similar  to  those  shown. 

At  this  point  it  is  possible  to  compare  the  calculated 
and  experimentally  observed  isotope  shifts  in  equilibrium 
F-F  distance,  I.R.  active  stretching  frequency  (v3) ,  and 
Raman  active  stretching  frequency  (v^) . 

A  small  negative  shift  upon  deuteration,  ARQ ,  has 
been  observed  in  the  equilibrium  F-F  distance  in  crystals 
of  sodium  and  potassium  bifluoride.  ^^2'  ARq  = 

-0.0024  A;20'29  NaHF2 /  ARQ  =  ■°*0046  A  20'28).  From  the 
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difference  between  the  minima  of  the  curves  Uq(R)  and 

H  o 

Uq (R) ,  ARq  has  been  estimated  as  -0.007  A,  which  is  in 

good  agreement  with  values  calculated  by  Alml0f2^ 

(-0.0058  A)  and  Jiang  and  Anderson^J  (-0.006  A). 

Using  the  harmonic  approximation,  and  have 

been  estimated: 

^  1516  cm"1  ;  1058  cm"1  ; 

giving:  Av3  -458  cm"1;  v^/v^  1.43.  Ketelaar  and 

95  96  -1  -l 

Vedder  '  have  measured  1450  cm  and  1023  cm  (Av3  = 

—  1  H  D 

-427  cm  ;  ^3/^3  =  1*42)  for  KHF3  and  KDF3  respectively. 

Finally,  the  Raman  active  transition  is  predicted 
to  have  a  small  positive  shift: 

H  677  cm  1  ;  680  cm  1  ;  Av^  ^  +  3  cm"1. 

The  calculation  of  Alml^f24  yielded  Av1  =  +6  cm”1. 

97  -1 

Dawson  has  experimentally  found  Av1  to  be  +1.25  cm 

(±0.25  cm  1)  for  the  potassium  salt. 

3.  Asymptotic  Values: 

Asymptotically,  each  of  the  two  wells  of  the  double 

minimum  potential  must  approach  the  potential  function 

for  a  free  HF  molecule.  The  fundamental  stretching 

—  1  98 

frequency  for  HF  is  3958.4  cm  ,  from  which  a  harmonic 
oscillator  force  constant,  k  =  0.56749  a.u.,  has  been 
calculated.  Because  the  adiabatic  representation  uses 
the  reduced  mass  for  HF2~>  m,  (2.6%  greater  than  the  HF 
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reduced  mass)  a  value  of  3908  cm”1  is  expected  for  the 
difference  between  the  asymptotic  values  of  ideal  curves 
U2  and  Uq.  At  large  R  values  (greater  than  7.0  a.u.) 
this  difference  actually  converges  to  about  3930  cm  1. 

The  asymptotic  curvature  of  the  potential  function  at 
each  minimum  is  clearly  a  little  high,  but  is  nevertheless 
a  good  representation. 

It  is  also  possible  to  estimate  asymptotic  values 

for  the  matrix  elements  {(i/h)P. .}. 

1  j 

At  any  R  point  where  the  potential  has  a  double 
minimum  (at  ±Zq(R))  it  may  be  possible  to  represent  the 
vibrational  Born-Oppenheimer  wavef unctions  by  a  series 
of  pairs  of  diatomic  molecular  vibrational  functions 
(Hermite  functions  in  the  harmonic  approximation)  centred 
on  +Zq  and  -zQ,  that  is,  a  two-centre  expansion: 

aB0(z;R)  =  l  a  .  (R)  (0  . "  (z )  +  (-1 )  U+ j  ,+ (z )  )  (4.27) 

4  4D  3  3 

where : 

*.+  (z)  =  exp (- (n (z±zQ) ) 2/Z) Hj (n (z±zQ)) 

e  W. (u)  (4.28a) 

3 

u  =  n(z±z0);  _CWjWkdu  =  6jk  (4.28b) 

Hj  (u)  is  a  normalized  Hermite  polynomial.  <J)j  is 
centred  on  “zq  ^nd  4)j+  is  centred  on  +zq* 

(4.27)  contains  the  term  (-1)P+:I 
(g,u)  symmetry  requirement  is  met: 


in  order  that  the 


« 
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ct®°(z;R)  =  (-1 )  y  a®°(-z;R)  (4.29) 

M  M 

n (R)  is,  in  general,  a  parameter  to  be  chosen  for 
convenience,  perhaps  defined  by  the  curvature  at  z  =  Zq. 

The  expansion  basis  is  non-orthogonal ,  and  as  the 
overlap  between  the  functions  {(j)^  ;  4>j  +  }  increases  (as 
Zq  -►  0,  the  barrier  height  goes  down)  it  becomes  increas¬ 
ingly  difficult  to  uniquely  define  coefficents  {a  .}. 

J 

At  large  R  distances,  however,  the  BO  states  become 
well  defined  degenerate  (g,u)  pairs  which  directly 
correspond  to  HF  molecular  functions.  If  the  HF  functions 
are  approximated  by  harmonic  states,  the  parameter  n  is 
now  fixed  using  the  force  constant,  k,  previously 
calculated : 

n  =  ( (km)1/2/h)1/2  (4.30a) 


which  gives 


X]  =  5.645  a  .u 


(dimension:  length  -'-) 


(4.30b) 


and  the  first  four  Born-Oppenheimer  states  are: 


a®°  =  (n/2) 1/2  ((j)0-  +  4>0+) 

a®0  =  (n/2) 1/2  (<t>0-  -  <t>0+) 
a®0  =  (n/2) 1/2  1*1  -  4>1+) 
a®0  =  (n/2) 1/2  +  <t>1+) 

(i/h) PQ2  (or  P13)  may  now  be  estimated, 
derivative  property  of  the  normalized  Hermite 


(4.31) 


Using  the 
functions , 


. 
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(4.28a)  gives: 


a/3R4» j-  =  ±nz0’  (R)  [  (j/2)1/2(J)  +  _1  -  (  (j  +  l)/2)1/2ci)++1] 


(i/fi)  PQ2  =  <a"u  la/BRla"  > 

/  (n/2)  (4>q+4>q)  HZq  [2  1^2  ( ^ o-1-^ 0 )  ”  ^2+(^  ]  dz 


(4.32) 


i.e.,  ( i/fi )  PQ2  ~  2  -n-z^(R) 


-1/2 


(4.33) 


However ,  zQ  (R)  ~  R/2;  Zq  (R)  ~  0.5 


and : 


(4.34) 


Figure  4.3  shows  that  (i/fi)P  2,  calculated  using 
numerically  generated  eigenfunctions  (solutions  to  (4.1)), 
converges  most  satisfactorily  to  the  asymptotic  value 
given  by  equation  (4.34). 

It  is  also  possible  to  estimate  asymptotic  values 

for  the  other  matrix  elements  in  this  way.  The  one-term 

harmonic  approximation  for  higher  states  will,  however, 

give  poorer  results;  for  example,  in  the  harmonic 

approximation  (i/fOP^  is  asymptotically  zero,  whereas 

figure  4.4  shows  it  to  be  about  0.2  a.u.  ^ .  This  merely 

BO 

indicates  the  asymptotic  anharmonicity  of  ,  which  may, 
however,  be  represented  by  a  series  of  the  harmonic 
functions  (as  in  equation  (4.27)),  the  coefficients  for 
which  will  be  asymptotically  constant,  leading  to  the 
constant  value  for  (i/fi)Pg^. 
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The  translational  origin  of  the  asymptotically 
non-vanishing  constant  value  for  (i/ft)PQ2  is  evident  in 
equation  (4.33).  In  this  derivation  the  proton  transla¬ 
tion  effect  has  appeared  as  the  term  z^ (R) ,  the  trans¬ 
lation  of  the  position  of  the  minimum  of  the  potential 
with  respect  to  the  centre  of  mass  of  the  heavy  nuclei. 

It  is  also  possible  to  obtain  an  estimate  of  the 
"correct  tailing"  of  the  modified  couplings  {C^}.  This 
"residual  coupling"  arises  from  a  distortion  of  the 
double  minimum  wells  as  the  bond  dissociates  (polariza¬ 
tion)  ,  and  may  be  estimated  by  replacing  the  operator 
[9V,/9R]  in  equation  (4.15)  by 

-L  Z 

[9V^/9R]g  for  z  <_  0 ,  and  by 

[ 9V-^/9R]  t  for  z  _>  0  (s  and  t  being  the  appropriate 
Jacobi  coordinates  z  ,z_  in  this  case:  see  figures  3.1 

A  jj 

and  2.2). 


The  resulting  matrix  elements  follow  the  curves 
{C . . (R) }  very  closely  for  R  >  5.4  a .  u .  ,  and  by  ~1 . 0  a . u . 

they  are  virtually  identical. 

This  shows  that  the  switching  function  corrections 
are  physically  reasonable:  they  remove  spurious  coupling 
due  only  to  translation,  while  retaining  coupling  due  to 
distortion  of  the  local  potential  wells. 


CHAPTER  FIVE 


SOLUTION  OF  THE  COUPLED  EQUATIONS 
A.  Eigenvalue  Condition. 

The  eigenvalue  problem  consists  of  the  search  for 
particular  values  of  E,  and  the  components  of  a  vector 
F (R) ,  which  satisfy  the  matrix  equation: 

{  (-h2/2y)  (d/dRI+C (R) ) 2  +  U(R)-EI}F(R)  =  0  (5.1) 

where  I_  is  the  unit  matrix,  U  is  a  diagonal  matrix  with 
elements  (Un(R)},  and  C  is  a  skew-symmetric  matrix  with 
elements  {C^  (R)  }  (figures  4.2  to  4.8). 

In  a  manner  similar  to  that  of  Chapter  Four  (section 
A) ,  left  and  right  asymptotic  regions  (LAR,  RAR)  are 
defined  to  be  regions  of  small  and  large  R  respectively, 
far  from  the  turning  points  of  the  curves  (Un(R)}. 

In  either  asymptotic  region,  C(R)  approaches  zero, 
and  the  nth  component  of  a  solution  vector  F(R)  (for 
arbitrary  E)  then  behaves  as  a  solution  to  the  uncoupled 
one-dimensional  Schrodinger  equation  with  potential  Un (R) . 
There  are,  therefore,  two  linearly  independent  solutions 
for  each  component  of  a  general  solution  vector  F (R) . 

This  may  also  be  expressed  in  the  following  way: 

there  exist  2N  linearly  independent  solution  vectors 

for  an  arbitrary  energy  E,  where  N  is  the  dimension 
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of  the  matrix  equation.  The  general  solution  is  a 
linear  combination  of  this  set. 
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For  an  eigenvalue  E  ^ ,  the  components  of  the  correspond¬ 
ing  vector  f j (R) ,  { f j R (R) ;  n  =  1,  2,  N},  must  each  be 

continuous  and  everywhere  'regular',  that  is  they  satisfy 
the  boundary  conditions: 


Lim{f.  (R)  }  =  0  as  R  -►  0  or  00 


(5.2) 


(5.2)  eliminates  N  of  the  independent  solutions 
(those  with  components  growing  exponentially  away  from 
the  turning  points) . 

For  any  E,  N  linearly  independent  solution  vectors 
may  be  constructed  which  are  regular  in  LAR,  and  N  may 
be  constructed  which  are  regular  in  RAR.  The  following 
set  may  be  used  at  an  initial  point  in  either  LAR  (R^) / 
or  RAR  (Rr) : 


1 — 

1 — 1 

1 _ 

t 

1 
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1 _ 
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1 _ 
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• 

• 
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• 

• 

• 

1 

0 

_ 1 

1 — 

0 

1 _ 

1— 1 

Clearly  this  set,  {F.k(R)?  k  -  1,  2,  .*.  N},  is 
linearly  independent  at  R^  or  R^.  Initial  values  for 
derivatives  may  be  fixed  by  the  JWKB  criterion  discussed 
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in  Chapter  Four,  although,  for  the  same  reasons,  the 
numerical  values  are  unimportant. 

The  solution  vectors , ( 5 . 3 ) ,  may  be  numerically 

propagated  from  an  asymptotic  region  towards  the  turning 

£ 

points  of  the  curves  (Un(R)}.  Let  (£n(R);  n  =  1,  2,  ...  N} 
and  n  =  1,  N}  be,  respectively,  left  and  right 

regular  linearly  independent  sets  of  solution  vectors. 

For  a  general  energy  E,  no  linear  combination  of  the  set 
{F  }  is  identical  to  any  linear  combination  of  the  set  {F'  }. 

It  is  only  when  E  is  an  eigenvalue  that  the  two  sets  may 
be  matched. 

Thus,  the  eigenvalue  condition  is: 

l  {dV  +  drFr}  =  0 
L ,  n— n  n— n  — 

n=l 


and 


l  {d  V  +  drFr'}  =  0 
L  ,  n— n  n— n  — 

n=l 


(5.4) 


(5.4)  implies  the  existence  of  the  2N  coefficients 

{d£,  dr;  n  =  1,  2,  . . .  N} ,  which  requires  the  vanishing 

n  n 

of  the  2N  X  2N  determinant: 
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Nl  *  *  * 
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Nl  *  *  * 

• 
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F 

NN 

The  principle  for  the  solution  of  the  eigenvalue 
problem  is,  therefore,  merely  a  generalization  of  the 
method  of  Chapter  Four.  The  energy  dependent  function 
W(E),  defined  by  (5.5),  may  be  called  a  generalized 
Wronskian  (for  N  =  1,  it  reduces  to  equation  (4.9)). 

The  problem  appears  to  be  trivial:  start  a  numerical 
integration  of  the  coupled  differential  equations,  (5.1), 
using  initialized  linearly  independent  column  vectors 
given  by  (5.3);  propagate  from  LAR  and  RAR  to  a  suitable 
matching  point,  in  a  classically  allowed  region;  calculate 
W(E)  given  by  (5.5),  and  numerically  determine  zeros.  The 
zeros  of  W(E)  are  eigenvalues,  { E ^ } . 

Two  main  areas  of  difficulty  arise:  linear  dependency 
and  nodes.  Before  discussing  them  it  is  necessary  to 
introduce  the  concept  of  open  and  closed  channels. 


4 
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B.  Open  and  Closed  Channels. 

1.  Zero  Coupling. 

If  C  is  everywhere  zero,  (5.1)  reduces  to  N  independent 
one-dimensional  Schroedinger  equations  for  the  motion  of 
the  reduced  mass  y  in  any  of  the  N  potential  curves 
(Un(R)}.  This  is  precisely  the  Stepanov  (or  Born- 
Oppenheimer)  approximation.  It  may  also  be  called  the 
one-channel  approximation,  as  wavef unctions  associated 
with  a  particular  potential  curve  evolve  independently 
of  functions  associated  with  all  other  curves. 

If  eigenvalues  are  sought  for  values  of  E  which  are 
less  than  the  minimum  of  the  curve  U-^  (R)  ,  there  can  only 
be  solutions  for  the  curve  UQ(R).  There  is  said  to  be 
one  OPEN  channel,  and  (N-l)  CLOSED  channels.  If,  however, 

E  is  greater  than  the  minimum  of  (R) ,  but  less  than  the 
minimum  of  U2 (R) ,  there  are  now  two  open  channels,  and 
(N-2)  closed  channels.  Eigenvalues  associated  with  U1 (R) 
and  UQ(R)  are  possible  in  this  case. 

In  this  simple,  zero  coupling,  case  it  is  easy  to 
distinguish  between  eigenvalues  associated  with  and 
,  simply  by  initializing  vectors  with  a  'signal'  (non¬ 
zero  value)  in  the  appropriate  channel;  for  example, 
of  equation  (5.3)  could  only  generate  eigenvalues  of 
U  (R) ,  and  F~  could  only  give  eigenvalues  of  U1 (R) . 

In  the  zero  coupling  case,  no  signal  can  ever  occur 
in  any  channel  which  is  initially  zero. 
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2.  Non-Zero  Coupling. 

Consider  a  two-channel  situation  (that  is,  the 
expansion  basis  is  limited  to  two  terms;  N=2)  with  non¬ 
zero  coupling.  The  upper  channel  is  again  said  to  be 
closed  if  E  is  everywhere  less  than  (R) .  In  order  to 
find  an  eigenvalue,  two  linearly  independent  solution 
vectors  must  be  generated  in  each  asymptotic  region,  and 
propagated  to  a  classically  allowed  region  of  UQ (R) . 
Although  a  vector  of  the  form  (equation  (5.3))  initially 
has  no  signal  in  the  closed  channel,  this  component  grows 
exponentially  during  the  numerical  integration  procedure. 

As  integration  continues  into  the  classically  allowed 
region  of  UQ(R),  the  closed  channel  signal  begins  to 
dominate  the  open  channel  signal.  The  result  is  that, 
despite  initializing  two  linearly  independent  vectors  of 
the  form  and  F 2 ,  they  rapidly  become  linearly  dependent 
as  integration  proceeds  through  a  region  of  non-zero 
coupling  in  which  at  least  one  channel  is  closed.  This 
problem  always  occurs  for  (N  >  1) -channel  coupling.  It 
has  been  found  that  two  initially  independent  vectors  may 
become  identical  in  their  components  to  at  least  eight 
significant  digits. 

The  existence  of  linear  dependence  among  either  the 

members  of  {F£},  or  of  {Fr},  invalidates  the  eigenvalue 

— n  — n 

condition  (5.4).  If  there  is  such  a  linear  dependency, 
the  determinant  W(E),  (5.5),  vanishes  for  all  E.  This 


. 
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tendency  to  linear  dependence  has  been  overcome  using  a 

reorthogonalization  procedure  (originated  by  Riley  and 
99 

Kuppermann;  also  described  and  applied  by  Wu  and 
Levine^O)  which  will  be  described  in  the  next  section. 

It  is  worth  noting  that  for  the  two  channel  case, 
with  both  channels  open,  it  may  become  difficult,  or 
impossible,  to  associate  any  eigenvalue  with  a  particular 
potential  surface.  For  a  weak  coupling  case,  however,  it 
is  expected  that  a  one-to-one  correspondence  of  eigen¬ 
values  and  potential  curves  will  be  maintained.  This 
will  be  discussed  in  greater  detail  with  reference  to 
the  results  obtained  for  the  HF2-  system. 

C.  Linear  Dependence. 

1.  Reorthogonalization. 

The  tendency  of  a  set  of  left  or  right  regular 

solution  vectors,  {F  ;  n  =  1,  2,  - N},  to  become  linearly 

— n 

dependent  may  be  compared  to  the  gradual  rotation  into 
one  plane  of  the  physical  vectors,  A,  B,  C,  initially 
chosen  to  be  orthogonal  (perpendicular  in  space) .  The 
three  vectors  may  be  used  to  define  a  parallelepiped, 
the  volume  of  which  is  given  by  the  triple  product: 
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A.BxC  = 


A 

x 


A 

y 


B 

X 


B 

y 


v 


(5.6) 


If  the  vectors  lie  in  a  plane  (V=0) ,  there  is  a  linear 
dependency  relationship  between  them.  A  necessary  condi¬ 
tion  for  linear  dependence,  therefore,  is  that  the 
determinant  V  shall  vanish.  If  the  three  vectors  are  not 
orthogonal,  but  do  not  lie  in  a  plane,  they  may  be  rotated 
back  to  a  position  of  mutual  perpendicularity.  This  is 
the  process  of  reorthogonalization :  a  particular  choice 
of  a  linear  combination  of  the  vectors. 

The  linear  independence  of  a  set  of  N  N-dimensional 
solution  vectors  may  be  maintained  by  periodic  reortho¬ 
gonalization  at  particular  points,  {R^},  in  the  numerical 
integration . 

The  NxN  determinant  V  (analogous  to  (5.6)),  is 

constructed,  with  each  column  being  one  of  the  vectors 

F  .  Each  column  is  'normalized'  by  dividing  the  com- 
— n 

ponents  of  that  column  by  the  square  root  of  the  sum  of 
the  squares.  The  value  of  this  determinant  is  initially 
one  (equation  (5.3)),  and  is  found  to  decrease  smoothly 
to  zero.  There  is  a  complication,  due  to  the  presence 
of  'nodes',  which  will  be  discussed  later.  If  the  value 
of  V  falls  below  0.95  (chosen  to  maintain  fairly  rigorous 


t 

' 
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orthogonality) ,  the  reorthogonalization  procedure  is 
instituted . 

Let  F(R^)  be  the  matrix  formed  by  the  columns  of 
the  vectors  (Fn(R^)}  at  a  reorthogonalization  point  R^ . 
Reorthogonalization  consists  of  finding  the  matrix  M  such 
that : 


F  (Rj_)  •M(Ri)  =  I  (5.7) 

where  M  is  the  inverse  of  F. 

At  the  point  R^ ,  new  derivative  vectors  must  also  be 
specified.  Although  the  order  of  the  matrix  multiplica¬ 
tion  in  (5.7)  is  unimportant,  it  represents  a  particular 
linear  combination  of  solutions  to  the  differential 
equation  (5.1),  and  the  correct  order  must  be  maintained 
in  redefining  the  derivatives: 

F ' (R.)*M(R.)  =  F'  (R. )  (5.8) 

—  l  —  l  —new  l 

2.  Extended  Definition. 

In  order  to  define  a  continuous  eigenfunction,  it 
is  necessary  to  apply  the  rotation  matrix  M(R^)  to  all 
points  in  the  numerical  integration  which  occurred  before 
the  reorthogonalization  point  R^.  That  is,  for  continuity, 
the  same  linear  combination  of  solution  vectors  must  be 
used  at  all  points. 

Suppose  the  matrices  {F^U^)}  have  been  stored  ( after 
reorthogonalization  where  applicable)  at  all  integration 
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points  {Rj};  this  set  has  discontinuities  at  the  reortho- 

gonalization  points  {R.  }.  A  new,  continuous,  set  {F  (R.)} 

i  —new  j 

is  constructed  in  the  following  way: 

Rnew(R)  =  £old(R)^(Ri>  (Ri-l^R<Ri>  (5.9) 

— new^Ri^  ^as  a^rea(^y  been  calculated  as  I_,  and  the 
point  R^_-^  used  a  different  matrix,  M(R^_^)  ,  to  produce 

^old  ^Ri-1  ^ =— '  so  the  defini'ti°n  of  F.new(R)  may  now  be 
extended : 


F 

—new 


(R) 


(Foid(R) •M(Ri_1) ) *M(Ri)  (Ri_2^R<Ri_1)  (5.10) 


and 


F 

—new 


(R) 


=  (Roid(R) *-(Ri-2} }  *-(Ri-l) '-(Ri} 


(R 


i_3<R<Ri_2) 


(5.11) 


In  this  way,  F(R)  is  redefined  (after  numerical 
integration  is  completed)  as  a  continuous  function  from 
the  matching  point  (where  a  final  reorthogonalization  is 
performed)  back  to  the  initial  point  in  either  asymptotic 
region. 


3.  Nodes. 

The  vanishing  of  the  determinant  V  is  a  necessary 
condition  for  the  linear  dependence  of  N  solution  vectors. 
It  is  not,  however,  a  sufficient  condition  to  ensure 
their  linear  dependence. 
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There  exist,  in  general,  particular  points,  {R  },  at 

P 

which  some  linear  combination  of  the  linearly  independent 

vector  solutions  to  equation  (5.1)  is  the  null  vector. 

These  points  are  called  nodes. 

The  matrices  {F(R  )}  are  then  singular,  but  the 

P 

derivative  matrices,  {F' (R  )},  are  not.  In  the  region 

P 

of  a  node,  therefore,  the  derivative  vectors  may  be 
reorthogonalized ,  if  necessary. 

There  are  also  points,  {R  },  at  which  the  determinants, 
{VP},  of  the  matrices  {F’ (R  )},  vanish.  These  points  may 
be  called  'derivative  nodes'. 

The  situation  may  be  summarized  in  the  following  way: 
in  general,  three  things  happen  as  the  numerical  integra¬ 
tion  progresses, 

(i)  initially  linearly  independent  solution  vectors 
become  dependent;  this  is  indicated  by  a  gradual 
decrease  in  both  |v|  and  | VP | ; 

(ii)  'nodes  of  V'  occur,  in  the  region  of  which  |v| 
diminishes  rapidly,  but  | VP |  remains  stable; 
derivative  vectors  may  be  reorthogonalized; 

(iii )  'nodes  of  VP'  occur;  | VP |  diminishes  rapidly, 
but  | V |  remains  stable;  vectors  may  be 
reorthogonalized . 

It  is  quite  often  unnecessary  to  reorthogonalize  in 
the  region  of  either  kind  of  node,  as  it  is  only  when  |v| 
and  I  VP  I  are  s'imu'ltaneousZy  small  that  the  solutions  are 
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truly  becoming  linearly  dependent. 

D.  Symmetric  Potential  Parity  Separation. 

If  V-^(z;R)  is  symmetric  in  z  space,  it  has  been 
shown  (Chapter  Four,  section  B5)  that  there  is  no  coupling 
between  Born-Oppenheimer  states  of  opposite  parity.  The 
coupled  equations  (5.1),  may  be  solved  independently  for 
each  parity,  in  this  case. 

Components  ^Fn j /  n  =  0,  2,  4  . .},  of  a  vector  F j , 
behave  independently  of  components  (Fr j ;  n  =  1,  3,  5  ..} 
(the  matrix  C (R)  is  block-factored  according  to  parity) , 
and  therefore  there  exists  a  set  of  eigenvalues  {E^}, 
associated  with  even  parity  curves  (Un(R);  n  =  0,  2,  ..}, 
and  a  completely  independent  set  of  eigenvalues  associated 
with  odd  parity  curves,  (Un(R);  n  =  1,  3,  ..}. 

E.  Eigenfunctions. 


1.  Construction: 

An  eigenvalue  Ej  has  an  associated  eigenfunction: 

N 

iK(z,R)  =  l  f  .  .  (R)a  (z;R)  .  (5.12) 

3  1=0  ^  * 


Let  the  vector  with  components  be  (R) • 

For  the  symmetric  potential  case,  the  eigenfunctions 
are  parity-separated,  and  (5.12)  becomes: 


(2N-2+P) 
lK(z,R)  =  l 
3  i=P,P+2  3X' 


(5.13) 
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where  P  =  0,  1  for  even,  odd  parity. 

It  is  necessary  to  distinguish  between  an  eigen¬ 
vector  3  f  ^ ,  and  a  set  of  N  linearly  independent  solution 

vectors,  {F,  ;  k  =  1,  2,  ..  N}.  The  eigenvector  f.  must 
— k  —j 

£ 

be  constructed  from  a  set  (F^(R)},  for  R  £  Rm,  (where  Rm 
is  the  matching  point,  at  which  W(Ej),  equation  (5.5), 

v- 

has  been  found  to  vanish)  and  a  set  (Fk(R)},  for  R  £  Rm. 
It  will  be  shown  that  reorthogonalization  of  {F^},  and 
{Fk},  at  Rm,  leads  to  a  simplification  in  the  eigenvector 
construction,  and  so  Rm  is  chosen  between  nodes  of  the 
determinant  V  (section  C3  of  this  chapter) .  Reortho¬ 
gonalization  at  the  matching  point  leads  to  a  matching 
of  the  magnitudes  of  left  and  right  solution  vectors. 

For  R  £  Rm,  the  eigenvector  is  a  particular  linear 
combination  of  the  set  {F^.},  but  for  R  £  Rm  it  is,  in 
general,  another  particular  linear  combination  of  the 


set  {F^}. 

That  is: 

f  j  (R) 

N  i  z 

=  I  dS(R) 

k=l  k 

v  | 

N 

(5.14) 

f  j  (R) 

=  l  d£f£(R) 
k=l  k  K 

3 

1  V 

In  order  to  define  the  (as  yet  unnormalized)  eigen¬ 
vector  f j (R) ,  for  all  R,  the  set  {d^;  d£}  must  be  known. 
At  Rm: 

f  .  (Rm)  =  l  d^F^ (Rm)  =  l  dk£k(Rm) 

3  k  k 


(5.15) 


' 
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and  the  reorthogonalization  gives: 


f  j  (Rm) 


rH 

+  d^ 

2 

- 1 

o 

_ 1 

+  • •  +  dN 

0 

0 

1 

0 

• 

• 

• 

• 

• 

• 

i - 

O 

» 

1 - 

o 

1 _ 

- 1 

r— H 

_  1 

'  l' 

+  dr 

'  o" 

+  . .  +  dr 

0 

2 

N 

0 

1 

0 

• 

• 

• 

• 

• 

• 

0 

0 

1 

(5.16) 


Thus  the  simplification: 


d 


l 

k 


(k  =  1,2, . .N) 


(5.17) 


The  coefficients  {d^.}  are  now  fixed  for  all  R  (providing 
appropriate  reorthogonalization  matrices  are  applied  in 
appropriate  regions  of  R) . 

Derivative  matching  provides  the  constraint  needed  to 
calculate  the  {d^}: 

l  dk(F^'  -  f£')  =  O  (5.18) 

k 

(5.18)  may  be  written  in  matrix  form: 


A .  d 


0 


(5.19) 


where : 


d 


1 


A  .  .  =  F .  .  -  F .  . 
1J  1] 


d 


d 


2 


d 


N 


(5.19)  is  an  NxN  homogeneous  system  in  the  N  unknowns, 
{d^},  and  can  only  yield  relative  values.  All  the  coef¬ 
ficients  may,  for  example,  be  calculated  relative  to  d^: 


(i  =  2,3, . .N) 


(5.20) 


Now  (5.19)  provides  N  ways  of  defining  the  (N-l) 

ratios  {a^}.  In  general,  for  large  matrices,  A_,  this  may 

become  a  sensitive  numerical  procedure,  but  for  the 

maximum  dimension  used  in  this  study  (N=3) ,  no  significant 

numerical  instability  has  been  found  (for  N=3,  all  three 

5 

ways  of  defining  a 2  and  a^  agree  to  1  part  m  10  ,  in  the 
worst  case,  over  the  whole  range  of  eigenvalues  calculated) . 

2.  Three  Channel  Example: 

The  three  channel  eigenfunctions  may  be  constructed 
using  (5.14) : 


f  .  (R)  =  dx  F^ttO  + 


a2  F12^  +  a3  F13 


F23  (R)  "  (ni1*1") 


« 


: 

* 

I 
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F11 <R) 

+  a2 

fI2(r) 

+  a3 

FI3(R) 

F21(R) 

F22  <R) 

F23  <R) 

>  (R>Rm) 

F31(R). 

F32(R). 

F33  (R). 

4 

(5.21) 

and  are  known  from  (5.19),  and  d^  is  to  be  fixed  by 
normalization . 


3.  Normalization: 

Given  that  the  components  {f.  },  of  the  eigenvector 

1  ~ 

f_j  (R)  ,  are  now  defined  continuously  throughout  R,  the 
constant  d^  is  determined  in  the  usual  way.  An  eigen¬ 
function  must  satisfy: 


.00  .00 


/  /  lj;  .  ( z  ,  R)  dzdR  =  1 


(5.22) 


—  00 


or ,  using  (5.12)  : 


/Co  r  ^  f  a  )  2dzdR  =  1  (5.23) 

0  -oo  £  ^  36 

Using  the  orthonormality  of  the  set  { ot ^ }  in  z  space, 
(5.23)  reduces  to  the  requirement: 


o  (£fj5,>dR  =  1 
£ 


(5.24) 


F.  IR  Relative  Intensities: 

The  integrated  absorptivity  associated  with  an  IR 
active  transition  from  the  ground  state  ijjg  ^  (the  super¬ 
script  indicates  even  parity,  P=0)  to  any,  odd  parity, 
0xcited  state  ijj  ^  ^  ,  is  shown  by  Wilson,  Decius  and  Cross 
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(reference  57,  p.  163)  to  be: 


/ Kdv  =  8 tt 3  v  (N*°*  -N  )  (Tn°)  2 
-  no  0  n 

3hc 


(5.25) 


(P) 


vnQ  is  the  frequency  for  the  transition,  N  '  the 


number  of  molecules  per  unit  volume  in  state  j ,  and  T 


no 


is  the  transition  moment  defined  by  (5.26)  (see,  for 
example,  Pauling  and  Wilson:  reference  62,  p.  299): 


Tno  E  f°°  j.™ 
o  _ 


(z,R)  y(z,R)  (z,R)dzdR; 


(5.26) 


y(z,R)  is  the  dipole  moment. 

Transitions  from  the  initial  state  ip 


( 0 ) 


( ' hot  bands ' ) 


will  not  be  considered  (the  occupancy  of  the  first  excited 
state  of  HF2~  at  90°K  is  about  10  ^  x  that  of  the  ground 


state) . 


The  intensity  of  the  transition  ip 


( 0 ) 
0 


(1) 


n 


,  relative 


to  the  intensity  of  the  ’  f undamentalT  transition  ip 
ipQ^  t  is  therefore  given  by: 


(0) 

0 


RI 


no  _ 


v. 


no 


v 


oo 


,no 


,oo 


(5.27) 


for  the  case  << 

n  0 


Using  (5.13)  the  transition  moment  integrals  become: 


,no 


r  /  (y (z ,R)  I 

0  “°°  K=  0 , 2  ,  . 


f0KaK 


I 

J=1 , 3 , . 


f  Ta,}dzdR 
n  J  J 


(5.28) 
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or 


no  r*  r  oo 

T  =  l  {/  f  OK  (R)  fnJ  (R)  [Zco  “K(z  ;R)  V  U  ’ R)  “j  (z;R)dz]dR} 
K'J  °  (5.29) 

A  Taylor  expansion  of  y(z,R)  about  the  point  z  =  0, 


.oo 


R  =  Rq  (the  position  of  the  minimum  of  Ug (R) )  is: 


y(z,R)  =  y (0 , Rg )  +  (R-Rg)3y/9R  +  z-8y/9z 

+  1/2 [ (R-Rg) 292y/9R2  +  2 (R-R  ) z92y/9R9z 
+  z  9  y/9z  ] 

+  higher  order  terms.  (5.30) 

(all  partial  derivatives  in  (5.30)  are  evaluated  at  z  =  0, 

R  =  Rg) . 

For  the  symmetric  case  the  dipole  moment  is  an  odd 
function  of  z(y(z,R)  =  -y(-z,R)),  and  all  coefficients  of 
even  powers  of  z  are  zero.  To  second  order,  the  effective 
dipole  moment  operator  may  thus  be  written: 


y  ^  A  (1  +  B  (R-Rg)  )  z 


(5.31) 


with : 


A  E  9y/9z  ;  A. B  E  92y/9R9z  . 


B |  is  expected  to  be  <1. 

A  second  order  estimate  of  the  relative  intensities 


is  then  given  by: 


. 
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no 


v 

no 


v 

oo 


;°{  £  f0KfnJIKJ}tl+B{R-R0,1'dR 

K ,  J 

I  f 0Kf 0JIKJ^  *-1+B  (R_R0^  *dR 
J\  /  J 


2 


(5.32) 


where : 


ikj(r) 


—  oo 


a(z;R)zaT(z;R)dz 
J\  J 


(5.33) 


For  three  channel  coupling,  there  are  nine  unique 
integrals  {I  }.  These  have  been  calculated  numerically 
(cf.  the  matrix  elements  { j }  of  Chapter  Four). 

The  effects  of  electrical  anharmonicity  may  be 
investigated  by  taking  |b|  >  0  in  equation  (5.32).  It 
will  be  shown  (Chapter  Six)  that,  at  least  for  HF2-, 
the  relative  intensities  of  combination  bands  of  the  funda¬ 
mental  (v^  +  n'v^')  and  overtone  ( ' 3v^  +  n'v^')/  are  not 
strongly  dependent  upon  the  value  of  the  parameter  B. 

That  is,  it  is  the  non-orthogonality  of  the  functions 

{fn„(R),  f  ,(R);  K  =  0,  2,  4,  ..  J  =  1,  3,  5,  ..}  (mechanical 
u  n  j 

anharmonicity  and  non-adiabatic  coupling  effects)  rather 
than  electrical  anharmonicity,  which  is  the  principal 
cause  of  the  broadness  and  structure  of  the  v3  band. 

G.  Curve  Fits  for  HF^-: 

The  sets  of  curves  (C^j(R)}/  (Un(R)}  and  (lKJ(R)} 
have  all  been  calculated  pointwise  in  R,  at  0.05  a.u. 
intervals,  from  3.80  a.u.  to  7.00  a.u.  It  is  necessary 
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to  have  continuous  representations  of  these  functions 
throughout  the  region  of  R  in  which  numerical  solutions 
to  (5.1)  are  generated. 

2 

Expansion  of  the  operator  (d/dR  +  C (R) )  ,  of 

equation  (5.1),  leads  to  a  term  involving  C*  (R) .  For 
this  reason,  it  is  convenient  to  represent  each  matrix 
element  by  a  single,  differentiable,  functional  form, 
throughout  R. 

Since,  however,  a  knowledge  of  the  derivatives  of 

{U  (R) } ,  and  {I  (R) } ,  is  not  explicitly  required,  they 
n  j\j 

have  been  represented  in  a  "piecewise"  manner;  first 
derivatives  are  matched  at  the  boundaries  of  each  region. 

1.  Matrix  Elements. 

For  three  channel  coupling,  it  can  be  seen  from 
figures  4.3  to  4.8  that  the  curves  { C ^ j (R) }  fall  into  two 
classes : 

(i)  adjacent  state  coupling. 

Curves  CQ2,  C24,  C13  and  C35  all  remain  positive, 

and  are  essentially  distorted  single  peaks.  The  Gaussian 

form,  Aexp (B (R-Rq) 2) ,  has  been  found  to  decay  too  rapidly 

at  large  | R-RQ | ,  but  a  hyperbolic  secant  form, 

Asech2 (B (R-Rq) ) ,  is  quite  satisfactory.  An  asymmetric 

curve,  C. . (R)  (j  =  i  +  2) ,  may  be  fitted  to  the  polynomial- 
lj 


modified  form: 


. 


sEi 
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cij  (R)  -  Y(R)  5  PL(R-R0)sech2(B(R-R0)) 


(5.34) 


where : 


pl(r-ro} 


L 

l 

1=0 


a£  ^R  V 


Z 


The  constants  Rq  and  B  are  initially  estimated  from 

the  position  of  the  maximum  of  C^j (R) ,  and  the  width  at 

half  the  height,  respectively.  For  a  given  order  L,  the 

coefficients  {a.},  of  the  polynomial  PT ,  are  then  obtained 

36  i-j 

by  a  least-squares  procedure,  that  is  by  minimizing  S 
defined  by: 


S  = 


M 

l  [C  .(R  )-Y(R.)] 
k=l  13  k  K 


(5.35) 


where  M  is  the  number  of  points. 

S  is  then  further  minimized  by  varying  B  (for  a 
fixed  Rq) ,  and  recalculating  the  polynomial.  Finally, 

S  is  minimized  by  varying  RQ ,  and  iterating  through  B 
and  P  .  In  principle,  higher  order  polynomials  will  also 
reduce  S,  but  in  practice,  of  course,  round-off  error 
limits  the  improvement.  Polynomials  of  order  eight  to 
ten  have  been  used. 


(ii)  Second  neighbour  coupling. 

The  shape  of  the  curves  Cq^  and  resembles  that 

of  the  derivative  of  a  single  peak,  and  the  curves  may 
kg  fitted  to  the  polynomial— modified  derivative  of  (5.34)  . 


. 

* 

. 
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cij (R)  ~  Y(R)  E  PL(R-RQ)sech2 (B(R-Rq) )tanh(B(R-RQ) ) 


(5.36) 


Again  the  parameters  B  and  Rq  are  initialized,  the 
least-squares  polynomial  obtained,  and  S  minimized  by 
varying  B  and  Rq . 

2.  (U_ }  and  {lTrT}: 

Let  qn (R)  represent  any  one  of  these  functions.  The 
(qn)  are  all  smoothly  varying,  non-oscillatory  functions, 
and  have  all  been  approximated  by  the  quadratic  spline-fit 
method . 

In  this  method,  the  region  R^  <_  R  £  R^ ,  in  which  the 
function  q (R)  is  known  pointwise  as  (q(R^);  R^  =  1,2,..M}, 
is  divided  into  sub-regions.  In  each  of  these  regions, 
q (R)  is  approximated  by  a  different  quadratic  polynomial, 
with  the  continuity  requirement  that  the  first  derivatives 
are  matched  at  the  boundaries  between  each  region. 

The  procedure  is  initiated  using  three  points,  in  a 
region  where  the  curvature  of  q(R)  is  small,  to  define  the 
coefficients  of  a  quadratic  for  region  1,  (R) : 


3  3 

Q,  (R)  =  I  a,  {  it  (R-R.)  } 
1  i=l  1  j=l  ] 


(5.37) 


whence 


a . 

l 


-1 


(5.38) 


. 
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Thus  q (R)  ~  Q1(R)  (R  <  R3)  (5.39) 

For  all  the  curves  (qn(R)},  the  first  three  points, 

R  =  3.80,  3.85,  4.00  a.u.,  have  been  used  to  define 
region  1. 

All  subsequent  regions  are  bounded  by  two  points. 

The  quadratic  applicable  to  any  region  is  defined  by  these 
two  points,  and  the  additional,  derivative  matching, 
constraint.  For  example: 

Q2(R3}  =  Q1(R3}  =  q(R3} 

Q2(R4)  =  q(R4) 

and  Q2 (R3 )  =  Q{(r3)  (5.40) 

This  procedure  provides  a  representation  of  q(R), 
which  passes  through  all  the  points,  R-^  to  R^,  and  has  a 
continuous  first  derivative. 


. 


CHAPTER  SIX 


RESULTS  AND  DISCUSSION  FOR  HF2- 
A.  Eigenvalues  and  Eigenfunctions. 

1.  Uncoupled  States: 

In  this  one-channel  case,  there  is  only  one  term  in 
the  expansion  (equation  (5.12))  of  the  eigenfunction 
ipj(z,R),  corresponding  to  an  eigenvalue  E  ^  .  Consequently, 
the  eigenstate  may  be  characterized  by  an  additional 
index,  n,  referring  to  the  potential  curve  U^ (R) ,  to 
which  it  belongs.  Thus: 

|  j  ,n>  5  ip  (z,R)  =  f  jn  (R)  an  (z  ;R)  (j  =  0,l,2,..)  (6.1) 

The  uncoupled  eigenvectors  {f^n(R)},  which  satisfy 
equation  (5.1),  are,  with  parity  separation  (P=0,1): 


fjp(R) 

•  f  = 

'  —  jP+2 

'  0 

•  f  = 

'  — j  P+4 

0 

0 

f.  (?) 
jP+2 

0 

0 

0 

f .  (R) 

jP+4 

• 

• 

• 

• 

• 

• 

Each  curve  U  (R)  has,  therefore,  associated  with  it 

n 

a  set  of  eigenvectors  and  eigenvalues  (_fjn(R)/*  Ejn'* 
j  =  0,1,2..},  which  have  no  dependence  on  any  other  curve 

Vr). 
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Figures  6.1  and  6.2  show,  for  HF^-f  the  series  of 


eigenvalues  calculated  for  even  and  odd  parity  curves  (Un). 

For  a  particular  value  of  n,  the  assignment  of  an 

eigenvalue  E^n  to  the  curve  Un (R)  is  unambiguous,  even 

in  the  case  of  the  near  degeneracy  of  and  E^^  (figure 

6.2).  E_ ^  is  clearly  associated  with  U^,  because  it  is 

/  0  \ 

generated  by  the  vector  u  ^  j  ,  whereas  E^-^  is  generated 


by 


f 71 (R) 
0 


03 


This  may  seem  painfully  obvious,  but  it  has  an  important 
conceptual  consequence:  the  one-channel  approximation 
asserts  that  the  observation  of  a  transition  |  0 , 0>->  |  j  ,n>  , 
allows  a  definite  statement  to  be  made  concerning  the 
change  in  the  vibrational  Born-Oppenheimer  state  of  the 
proton  (actually  the  reduced  mass  m) .  With  final  state 
| 0 , 1 > ,  the  proton  is  said  to  be  in  its  first  excited  state; 
a  'fundamental'  transition,  with  frequency  v^,  correspond¬ 
ing  to  an  excitation  of  the  ' H-F  stretching  mode',  is  said 
to  have  occurred.  [Actually,  of  course,  the  change  in 
state  of  the  system  is  not  quite  the  same  as  a  simple 
excitation  of  the  proton,  because  the  function  f Q0 (R)  is 
not  the  same  as  fQ1 (R) :  they  are  both  ground  state  wave- 
functions  for  the  reduced  mass  y,  but  in  different 
potential  wells  ].  Similarly  the  transition  |0,0>->|0,3> 
is  referred  to  as  an  'overtone'  of  the  H— F  stretching 
frequency: 


. 
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FIGURE  6.1 

EVEN  PARITY  UNCOUPLED  EIGENSTATES 


0 
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R(a.u.) 
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FIGURE  6.2 

ODD  PARITY  UNCOUPLED  EIGENSTATES 
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20  - 
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In  this  way,  the  line  due  to  the  transition  |0,0>+ 
|j,l>,  may,  with  some  validity,  be  referred  to  as  a 
'combination  band',  of  the  form:  +  j'v^'.  Again,  the 

frequency  has  been  written  within  quotation  marks 
because  the  spacing  of  levels  (corresponding  to  energy 
states  of  the  reduced  mass  y,  in  the  one-channel  approx¬ 
imation)  in  the  curve  U-^  (R)  is  not  the  same  as  the  Raman 
active  transition  frequency  v-,  :  this  corresponds  to  the 
transition  |0,0>+  |l,0>  (excitation  of  the  ' F-F  stretching 
mode ' ) . 

To  summarize,  in  the  one-channel  approximation,  the 
terms:  'state  of  the  proton',  and  'combination  bands  of 

the  fundamental  (and  overtone  'Sv^')  with  v-^',  have 
definite  meaning. 

2.  Coupling  Effects: 

For  N  (>1)  -channels ,  the  total  wavefunction  ipj(z,R) 
cannot  be  said  to  correspond  to  any  definite  vibrational 
state  of  either  the  proton,  or  the  reduced  mass  y.  It  is 
not  now  strictly  possible  to  characterize  the  state 
by  an  index  n,  referring  to  a  particular  curve  Un (R) , 
since,  for  two  channels: 

iK  (P)  =  f  jp  (R)  ap  (z  ;R)  +  f  jp+2  (R)ap+2  (z;R)  .  (6.3) 

(the  parity  superscript,  P,  is  only  possible  for  a 
symmetric,  AHA,  system).  Eigenvectors  satisfying  (5.1)  are: 


£ 


-  . 
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fjp(R) 

f jP+2 


for  example: 


II 

o 

< 

<4-1 

l 

Hi 

O 

O 

II 

I — 1 

M-J 

•« 

foi 

-u 

f  02 

o 

ll 

* 

f  03 

(6.4) 


If  the  coupling  between  Born-Oppenheimer  states  P 

and  P+2  were  strong,  the  contributions  of  fjp(R)  and 

(P ) 

f  jP+2  (R)  '  to  t^ie  state  vector  ^  (z,R),  could  be 

comparable,  and  it  would  not  be  possible  to  describe  a 


transition  from  ip 


0 


^  to  ijjj  as  a  combination  band, 


V3  +  ov  as  an  overtone  '3v3‘  +  j'vp'.  It 


would  be 

meaningless  to  talk  of  the  'vibrational  state  of  the 


proton'  or  'the  vibrational  state  of  the  reduced  mass  y ' . 

The  vibrational  state  of  a  closely  coupled  system 

would,  therefore,  be  completely  characterized  by  the 

(P) 

function  ih  .  ,  without  reference  to  the  isolated  states 

3 

of  excitation  of  the  symmetric  (F-F) ,  and  asymmetric 

(H-F) ,  stretching  modes.  It  would  only  be  possible  to  describe 

(P)  (P) 

a  transition  as  symmetric  ( ip  .  ip,  )  ,  or  asymmetric 

J  K 

(ip  .  -*  ip  ^),  but  it  would  not  be  possible  to  ascribe 

J  k 

to  that  transition  an  excitation  of  only  the  H-F,  or  only 
the  F-F,  mode.  In  this  case,  the  idea  of  isolated 
stretching  modes  becomes  invalid,  and  the  motions  are 
truly  inseparable. 


rv 


■ 


. 
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Fortunately,  the  non-adiabatic  coupling  turns  out  to 
be  weak,  and  a  one-to-one  mapping  of  the  uncoupled  to 
coupled  eigenvalues  is  possible.  The  double  index  nota¬ 
tion,  |j,n>,  remains  valid,  with  the  understanding  that 
the  proton  is  'almost'  purely  in  the  nth  Born-Oppenheimer 
state . 

For  weak  coupling,  the  index  n  defines  the  dominant 
channel.  For  example,  with  two  channels,  either: 


/”  (f.D(R))2dR  >>  /  (f  ,D,,(R))2dR  =>  |j,P>  (6.5) 

o  J  r  o  J  ^ 

indicates  the  lower  channel  dominates,  or  the  reverse  of 
the  inequality  indicates  the  upper  channel  dominates: 


/  (f.D,,  (R))2dR  >>  /“  (f  ,C(R)  )2dR  =>  |  j  ,P+2>  (6.6) 

o  j  a  o  J  r 

If  (6.5)  holds,  the  wavefunction  may  be  characterized 
by  the  double  index  |j,P>,  and  the  eigenvalue  by  EjP* 

This  is  always  the  case,  for  HF2~,  when  the  upper  channel 
is  closed  (states  EQ0  to  E4Q,  figure  6.1;  EQ1  to  E^, 
figure  6.2). 


A  more  interesting  situation  occurs  when  both  channels 

(P) 

are  open.  Upon  finding  an  energy  E  =  E ^  ,  for  which 

the  generalized  Wronskian  (W(E)  (equation  (5.5)),  is  zero, 
it  is  then  only  by  an  examination  of  the  components  of 
the  eigenvector  f ^  ^  (R) ,  that  an  appropriate  second  index 
may  be  assigned,  using  the  criteria  (6.5)  and  (6.6). 


* 
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In  particular,  the  states  |7,1>  and  |0,3>,  which 
are  very  nearly  degenerate  in  the  one-channel  approxima¬ 
tion,  are  noticeably  separated  by  the  introduction  of 
coupling : 


N  Channels: 
1 
2 
3 


-E71> 


-1 

cm 


-2.2 


+  37.2 

+55.5  (6.7) 


Figure  (6.3)  shows,  for  two  channels,  the  behaviour 
of  W (E)  in  the  region  of  these  closely  spaced  eigenvalues. 

Even  and  odd  parity  eigenvalues  (E^p,  Ejp+2^'  0^)tainec^ 
from  calculations  using  one  to  three  channels,  are  shown 
in  tables  one  and  two  respectively. 

The  effect  of  non-adiabatic  coupling  is  clearly  quite 
small  for  all  cases  in  which  the  upper  channel  is  closed 
(less  than  1/2%) ,  and  as  expected,  the  higher  eigenvalues 
(upper  channel  open)  are  more  noticeably  affected  (for 
example,  and  eq3  change  by  about  1%)  . 

The  convergence  of  the  eigenvalues  is  excellent, 
particularly  for  those  associated  with  the  lower  channel 
of  each  parity. 

The  weaker  convergence  of  the  eigenvalues  (eq^  and 
E13)  associated  with  the  upper  surface  (R) ,  is  expected: 
for  such  states  the  three-channel  approximation  links 
only  one  higher  surface  (U^ (R) ) .  It  would  require  a 


v 
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FIGURE  6.3 

W(E)  x  io2  VS.  Ex  10  (A.U.) 

2  CHANNELS  :  E  .  E 

71  03 
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TABLE  ONE 


EVEN  PARITY  EIGENVALUES  (HF2~) 
(ATOMIC  UNITS) 


E.n  x  102 

j0 

j 

1  Chan . 

2  Chan. 

3  Chan. 

0 

0.524433 

0.526146 

0.526148 

1 

0.832564 

0.834527 

0.834538 

2 

1.13347 

1.13555 

1.13560 

3 

1.42397 

1.42511 

1.42531 

4 

1.69966 

1.69672 

5 

1.95641 

1.94562 

6 

2.19363 

2.17827 

E.«  x  102 

3  z 

j 

0 

2.01665 

2.03001 

1 


2.30661 


2.32513 


X 
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TABLE  TWO 

ODD  PARITY  EIGENVALUES  (HF2“) 
(ATOMIC  UNITS) 


E  . ,  x 

3 1 

102 

j 

1  Chan . 

2  Chan . 

3  Chan . 

0 

1.21456 

1.21878 

1.21879 

1 

1.48124 

1.48488 

1.48492 

2 

1.73719 

1.73937 

1.73952 

3 

1.98096 

1.98057 

1.98089 

4 

2.21185 

2.20857 

2.20895 

5 

2.43015 

2.42493 

2.42514 

6 

2.63689 

2.63055 

2.63061 

7 

*2.83334 

2.82596 

2.82601 

8 

3.02071 

3.01274 

3.01268 

E  .  x 
33 

102 

j 


0 

*2.83233 

2.84289 

2.85131 

1 

3.09917 

3.11243 

3.11734 

* 

Near  degeneracy:  E71;  E03  * 


'*'*  N 

' 

' 
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four- channel  calculation  to  match  the  convergence  of  the 
lower  surface  eigenvalues.  Considering  this  convergence, 
however,  it  can  be  expected  with  confidence  that  a  four- 
channel  calculation  would  not  significantly  alter  the 
three-channel  upper  eigenvalues  shown. 

Three-channel  calculations  have  not  been  performed 
for  the  higher  states  of  even  parity  (Table  one)  for  two 
reasons:  (i)  they  are  not  involved  in  the  IR-active 
vibrational  transitions;  (ii)  the  effects  of  three-channel 
coupling  can  be  inferred  from  the  odd  parity  set  (Table 
two)  . 

Figures  6.4  to  6.15  show  the  components  (from  calcula¬ 
tions  using  one  to  three  channels)  of  eigenvectors 

{f.  (R)  } ,  associated  with  wavef unctions  { .  (z,R)}.  Each 

— j  n  j  n 

figure  has  six  numbered  graphs:  1 1'  shows  the  single 

component  f^n(R)  (either  f^p  or  fjP+2'  depending  on  the 

curve  U  (R)  to  which  it  belongs)  obtained  from  a  one- 
n 

channel  calculation;  '2',  and  '3'  show  components  f^p(R), 

and  f._in(R)  from  a  two-channel  calculation;  '4',  '5', 

jP+2 

and  '6'  show  components  f^p(R),  f jP+2 (R) ,  and  f jp+4 
from  a  three-channel  calculation* 

These  pictures  clearly  illustrate  the  concept  of  'weak' 
coupling  (compare  the  relative  magnitudes  of  the  components 
of  an  (N>1) -channel  eigenvector),  and  the  consequent 
assignment  of  the  second  index,  n,  indicating  the  dominant 
channel,  or  as  an  alternative  viewpoint:  the  uncoupled 
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FIGURE  6.4 
PSI  00 
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FIGURE  6.5 
PSI  01 
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FIGURE  6.6 
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FIGURE  6.7 
PSI  21 
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FIGURE  6.8 
PSI  31 
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FIGURE  6.9 
PSI  41 
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FIGURE  6.10 
PSI  51 
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FIGURE  6.11 
PSI  61 
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FIGURE  6.12 
PSI  71 
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FIGURE  6.13 
PSI  03 
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FIGURE  6.14 
PSI  81 
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FIGURE  6.15 
PSI  13 
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eigenstate  |j,n>,  from  which  the  coupled  state  is  'derived' 
(compare  the  functional  forms  of  the  dominant  and  uncoupled 
components) . 

It  is  worth  discussing  these  wavefunction  pictures 
in  some  detail. 

Figure  6.4,  labelled  'PSI  00',  shows  the  components 
of  the  even  parity  ground  state  eigenvector  f_^^  (R)  . 

Graphs  1,  2,  and  4  show  the  open  channel  function  f qq (R) r 
in  the  one-,  two-,  and  three-channel  cases  respectively. 
This  function  is  evidently  very  similar  to  that  of  the 
ground  state  of  an  harmonic  oscillator,  and  it  is  signifi¬ 
cant,  as  an  indication  of  the  weakness  of  the  coupling, 
that  the  introduction  of  coupling  does  not  visibly  alter 
this  functional  form.  It  has  a  very  much  greater  maximum 
value  ( ~ 2 . 0 )  than  the  closed  channel  functions  f ^ (R)  r 
graphs  3  and  5  (maximum  -0.03),  and  f ^ (R) ,  graph  6 
(maximum  -0.0012).  In  this  case,  there  can  be  little 
doubt  that  the  state  may  be  labelled  |0,0>,  and  that  the 
proton  may  be  considered  to  be  'almost  totally'  in  the 
Born-Oppenheimer  ground  state.  Graphs  3  and  5,  for 
f Q2 (R)  (two-  and  three-channel  cases)  show  a  'dip'  as 
fQ0(R)  begins  to  grow,  after  which  f Q2 (R)  'follows'  the 
form  of  the  open  channel  function,  but  with  a  slight 
'delay'.  This  'dip'  followed  by  'delayed  mimicry'  appears 
to  be  characteristic  of  the  first  dosed,  channel  for  all 
subsequent  states:  |0,1>  to  |6,1>;  the  oscillatory 
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behaviour  of  the  closed  channel  function  is  essentially 
determined  by  that  of  the  open  channel  function  (a 
"driven  oscillator"  response) . 

Figure  6.5  justifies  the  odd  parity  ground  state 
labelling  |0,1>.  It  is  very  similar  to  the  diagram  for 
|0,0>,  and  requires  no  further  comment. 

The  states  |l,l>  to  |6,1>,  figures  6.6  to  6.11,  have 
in  common  that  the  upper  channels  are  closed .  The  lower 
channel  function  fj^(R),  dominates  in  all  cases,  and  has 
the  form  of  the  jth  state  of  a  perturbed  harmonic  oscil¬ 
lator.  Note  the  characteristic  'mimicry'  of  the  upper 
channels . 

Figures  6.12  and  6.13  show  the  first  cases  in  which 
the  first  upper  channel  is  open.  These  two  also  correspond 
to  the  case  of  near  degeneracy  in  the  one-channel  calcula¬ 
tion:  | 7 , 1 >  and  | 0 , 3 > .  Consider  the  two-channel 

situation.  The  functional  form  of  the  dominant  channel 
(figure  6.12,  graph  2:  lower  channel  for  | 7 , 1 > ;  figure 
6.13,  graph  3:  upper  channel  for  |0,3>)  clearly  indicates 
the  appropriate  pair  of  indices.  Once  again,  these 
functional  forms  are  not  visibly  altered  by  the  introduc¬ 
tion  of  two-  or  three-channel  coupling.  There  is,  however, 
a  strikingly  different  behaviour  in  the,  now  open,  weaker 
channel  (graph  3:  upper  channel  for  |7,1>;  graph  2: 
lower  channel  for  |0,3>).  The  dominated  function  now 
exhibits  an  independent  form,  rather  than  merely  mimicking 
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the  behaviour  of  the  dominant  channel  (as  in  the  case  of 
a  closed  upper  channel).  For  |7,1>,  the  lower  channel 
imposes  some  oscillation  upon  the  ground  state  form  of 
the  upper  channel,  while  for  |0,3>  the  upper  channel 
significantly  perturbs  the  seventh  harmonic  oscillator 
functional  form  of  the  lower  channel,  but  the  character¬ 
istics  of  the  weaker  channel  are  maintained.  In  the 
three-channel  calculations  (graphs  4  to  6  of  both  figures) 
the  weaker  open  channel  function  is  visibly  perturbed  by 
the  coupling  to  a  higher  channel.  The  closed  (third) 
channel  again  essentially  mimics  the  behaviour  of  the 
channel  immediately  below  it.  The  significant  'mixing' 
of  the  open  channel  functions  for  the  closely  degenerate 

.  ]_4 

states  |7,1>  and  |0,3>,  leads  to  a  'Fermi  resonance' 
effect,  which  will  be  discussed  in  greater  detail  in 
section  B  of  this  chapter. 

Figures  6.14  and  6.15,  for  |8,1>  and  |l,3>,  also  have 
two  open  channels.  One  is  clearly  dominant,  allowing  the 
distinct  labelling.  The  weaker  open  channel  again  evolves 
independently,  but  with  superimposition  of  the  dominant 
channel  form,  and  the  closed  channel  mimics  its  immediately 
inferior  neighbour. 
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3.  Summary. 

Calculations  using  one  to  three  channels  have  been 
shown  to  produce  convergent  eigenvalues.  This  in  itself 
implies  that  the  coupling  is  weak,  a  conclusion  which  has 
been  confirmed  by  a  discussion  of  the  relative  magnitudes 
and  functional  forms  of  the  components  of  the  eigenvectors. 
The  latter  considerations  have  also  justified  retention 
of  a  double  index  notation  for  the  eigenstates  in  the 
(N>1 ) -channel  situation.  Hence  the  terms  'excitation  of 
H-F  and  F-F  modes'  have  been  shown  to  adequately  (but 
approximately)  describe  transitions  between  these  eigen¬ 
states.  A  description  of  the  IR  spectrum  of  HF2~,  using 
the  concepts  of  combination  bands  and  overtones  of 
and  is  therefore  a  reasonable  physical  approximation. 


B.  Transition  Frequencies  and  Relative  Intensities. 

Table  three  shows  a  comparison  of  transition  frequencies 

(E.  -E,  » ,  cm”'*')  obtained  from  this  calculation,  using  one 
jn  kO 

95-97  101-114 

to  three  channels,  with  those  observed,  '  and 

24 

calculated  by  Alml^f. 

The  frequency  (v-^)  calculated  for  the  Raman  active 
transition,  |0,0>  ->  |l,0>,  is  about  10%  higher  than 
experimental  values,  but  is  within  3%  of  the  value  obtained 
by  Alml0f.  This  is  probably  a  reflection  of  the  fact 
that  the  ab  initio  calculation  tends  to  overestimate  the 
curvature  of  the  potential  surface  near  the  equilibrium 
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TABLE  THREE 

TRANSITION  FREQUENCIES  (cm”1) 


FROM  INITIAL  STATE  PSI  0 

Final  Assign.  # 

State  ^ 

,0: 

Channels 

2  3 

Observed3 

Alml^f ^ 

1,0 

vi 

676.3 

676.8 

676.8 

600,631 

660 

0,1 

v3 

1514.6 

1520.2 

1520.2 

1425-1550 

1497 

1,1 

V3+Vl 

2099.9 

2104.2 

2104.3 

2045-2060 

2070 

2,1 

V2V1 

2661.7 

2662.7 

2663.0 

2610-2670 

2632 

3,1 

V3V1 

3196.7 

3192.1 

3192.8 

3166-3240 

3226 

4,1 

V3+4V1 

3703.4 

3692.5 

3693.3 

3916? 

5,1 

v3+5v1 

4182.6 

4167.3 

4167.8 

•p 

6,1 

v3+6v1 

4636.3 

4618.6 

4618.7 

4750? 

7,1 

v3+7v1 

5067.5 

5047.5 

5047.6 

0,3 

3v3 

5065.2 

5084.7 

5103.1 

5085-5099 

5157 

rH 

00 

v3+8v1 

5478.7 

5457.4 

5457.3 

1,3 

3v3+V1 

5650.9 

5676.2 

5687.0 

5590-5650 

FROM  INITIAL  STATE  PSI  1 

,0: 

0,1 

V3"V1 

838.4 

843.3 

843.3 

880,855 

837 

aReferences  95-97,  101-114. 
^Reference  24. 


N 

• 
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7  24 

value,  '  and  also  considers  an  isolated  ion  with  no 
crystal  field  effect. 

Frequencies  calculated  for  the  lower  combination 
bands  (v^  +  j'v  j  =  0,..3)  agree  well  with  both  observed 
lines,  and  the  values  of  Alml^f.  The  higher  combination 
bands  (j  =  4.. 6)  do  not,  however,  correspond  in  a  1:1 
way  with  those  reported  by  Alml^f,  a  fact  which  is  not 
surprising  considering  the  behaviour  of  his  potential 
function  (Chapter  three) . 

Tables  four  to  six  show  relative  intensities  calculated 
using  equations  (5.32)  and  (5.33),  and  incorporating  the 
'electrical  anharmonicity  parameter'  B  (defined  in  equa¬ 
tions  (5.31)).  It  is  evident  from  these  tables  that 
electrical  anharmonicity  is  not  a  primary  cause  for  the 
relatively  high  intensities  of  the  combination  bands  of 
the  fundamental  (v^  +  j'v^')  and  overtone  ('3v^'  +  j'v^'). 
In  general,  a  negative  B  increases  the  relative  intensities 
of  combination  bands  while  decreasing  the  overtone 
intensity,  and  a  positive  B  decreases  combination  band 
intensity  while  increasing  the  intensity  of  the  overtone. 
The  unusually  large  relative  intensities  are  due  to  a 
vibrational  Franck-Condon  effect,  that  is,  the  overlap  of 
the  functions  (fgK(R),  fnJ(R)}. 

For  a  consideration  of  the  effects  on  the  intensities 
due  to  coupling,  the  case  in  which  B  is  zero  (Table  four) 
is  quite  representative.  The  relative  intensity  of  a 


-♦ 
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TABLE  FOUR 


RELATIVE  INTENSITIES: 

B  =  ZERO 

Final 

#  Channels 

State 

1 

2 

3 

rH 

K, 

o 

0 . 1000E+1 

0 . 1000E+1 

0 . 1000E+1 

rH 

i — 1 

0. 2849E+0 

0. 3148E+0 

0. 3148E+0 

2,1 

0 . 5514E-1 

0.5979E-1 

0 . 5976E-1 

3,1 

0.8067E-2 

0 . 8636E-2 

0 . 8623E-2 

rH 

0 . 9335E-3 

0 . 1089E-2 

0.1084E-2 

• — i 

in 

0 . 1036E-3 

0 . 1336E-3 

0 . 1347E-3 

6,1 

0.1375E-4 

0 . 1496E-4 

0 . 1529E-4 

7,1 

0.1792E-5 

0 . 2754E-4 

0 . 9730E-5 

0,3 

0. 6026E-2 

0 . 3666E-2 

0. 3857E-2 

i — 1 

00 

0 . 1136E-6 

0 . 3194E-5 

0 . 2399E-5 

1,3 

0. 3545E-2 

0. 3468E-2 

0 . 3571E-2 
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TABLE  FIVE 


RELATIVE  INTENSITIES:  B  NEGATIVE 


B  =  -0.2 

Final  #  Channels 


State 

1 

2 

3 

0,1 

0 . 1000E+1 

0 . 1000E+1 

0 . 1000E+1 

1,1 

0 . 3101E+0 

0 . 3405E+0 

0 . 3404E+0 

2,1 

0 . 6388E-1 

0 . 6920E-1 

0 . 6916E-1 

pH 

CO 

0 . 9979E-2 

0 . 1071E-1 

0 . 1069E-1 

1 — 1 

0.1253E-2 

0.1446E-2 

0.1440E-2 

1 — 1 

in 

0.1497E-3 

0 . 1908E-3 

0 . 1897E-3 

6,1 

0.2036E-4 

0 . 2268E-4 

0 . 2308E-4 

i — 1 

0.2746E-5 

0 . 2333E-4 

0.7377E-5 

0,3 

0 . 5723E-2 

0 . 3611E-2 

0 . 3796E-2 

8,1 

0 . 2235E-6 

0 . 2764E-5 

0 . 2020E-5 

1,3 

0. 3568E-2 

0 . 3565E-2 

0.367  OE-2 

B  =  -0.5 


i — 1 

o 

0 . 1000E+1 

0 . 1000E+1 

0 . 1000E+1 

1,1 

0. 3515E+0 

0. 3824E+0 

0. 3824E+0 

2,1 

0 . 7 879E-1 

0 . 8  519E-1 

0 . 8  515E-1 

rH 

K 

CO 

0 . 1337E-1 

0 . 1437E-1 

0.1435E-1 

rH 

0.1845E-2 

0 . 2104E-2 

0 . 2095E-2 

5,1 

0 . 2384E-3 

0 . 2954E-3 

0 . 2939E-3 

6,1 

0 . 3326E-4 

0 . 3789E-4 

0. 3836E-4 

rH 

r- 

0.4638E-5 

0.1746E-4 

0 . 4362E-5 

0,3 

0. 5266E-2 

0 . 3527E-2 

0.3702E-2 

00 

I-1 

0 . 4681E-6 

0 . 2156E-5 

0 . 1495E-5 

1,3 

0. 3604E-2 

0 . 3719E-2 

0.3827E-2 

!  1  L  .  *  * 
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TABLE  SIX 


RELATIVE  INTENSITIES:  B  POSITIVE 


B  =  +0.2 


Final 

State 

1 

#  Channels 

2 

3 

0,1 

0 . 1000E+1 

0 . 1000E+1 

0 . 1000E+1 

1,1 

0 . 2615E+0 

0 . 2908E+0 

0 . 2908E+0 

2,1 

0.4728E-1 

0 . 5131E-1 

0 . 5128E-1 

3,1 

0.6405E-2 

0 . 6837E-2 

0. 6826E-2 

4,1 

0 . 6684E-3 

0.7898E-3 

0 . 7860E-3 

5,1 

0 . 6700E-4 

0 . 9095E-4 

0 . 9030E-4 

6,1 

0 . 8565E-5 

0 . 8986E-5 

0.9257E-5 

7,1 

0 . 1060E-5 

0 . 3197E-4 

0 . 1233E-4 

0,3 

0.6328E-2 

0 . 3719E-2 

0 . 3917E-2 

8,1 

0 . 4220E-7 

0 . 3642E-5 

0.2798E-5 

1,3 

0 . 3523E-2 

0 . 3374E-2 

0.3477E-2 

+  0.5 

i — 1 

o 

0 . 1000E+1 

0 . 1000E+1 

0 . 1000E+1 

1,1 

0 . 2294E+0 

0 . 2577E+0 

0 . 2  577E+0 

2,1 

0 . 3697E-1 

0 . 4  017E-1 

0 . 4  015E-1 

3,1 

0 . 4342E-2 

0.4604E-2 

0 . 4596E-2 

4,1 

0. 3633E-3 

0 . 4418E-3 

0 . 4393E-3 

5,1 

0 . 2811E-4 

0 . 4205E-4 

0 . 4166E-4 

6,1 

0.3232E-5 

0 . 3022E-5 

0.3188E-5 

7,1 

0.3401E-6 

0 . 3897E-4 

0 . 16  63E-4 

0,3 

0. 6778E-2 

0.3798E-2 

0.4004E-2 

8,1 

0 . 2300E-9 

0.4343E-5 

0 . 3431E-5 

1,3 

0.3491E-2 

0 . 3241E-2 

0.3342E-2 

' 
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combination  band  +  j'v^'  (j  =  1  to  4),  is  increased  by 
about  10%  in  going  from  one  to  two  channels  (compared  to 
a  change  of  less  than  1%  in  the  eigenvalues) ,  and  is 
virtually  unaltered  (j  =  1  to  5)  by  the  expansion  to  three 
channels.  This  behaviour  is  a  good  indication  of  the 
convergence  of  the  wavef unctions  for  these  lower  states. 

There  is  a  strong  'Fermi  resonance'  effect  for  the 
closely  degenerate  pair  |7,1>  and  |0,3>:  the  combination 
band  +  7'v^'  shows  an  approximately  five-fold  increase 
in  intensity,  accompanied  by  about  a  30%  decrease  in 
intensity  for  the  overtone  This  degeneracy  is 

'accidental':  it  is  a  consequence  of  the  particular  form 

used  as  a  model  for  the  true  potential  function,  and  in 
the  following  discussion  of  the  IR  spectrum  of  KHF^,  it 
will  be  argued  that  the  degeneracy  is  not  in  fact  'real'. 
The  existence  of  such  Fermi  resonances  in  the  spectrum  of 
a  hydrogen  bonded  system  is,  however,  highly  probable, 
even  if  so  far  unobserved. 

The  IR  transmission  spectrum  of  KHF2  at  90°K,  obtained 
by  Cote  and  Thompson^01  (C&T  henceforth)  is  shown  in 
figure  6.16. 

The  upper  curve  (1100  to  2700  cm  ^)  shows  vertical 
lines  placed  at  the  maxima  of  the  bands  assigned  by  C&T 
to  v  3  +  j  *  v 1 *  ,  for  j  =  0,1,2.  The  height  of  each  line 
is  proportional  to  the  calculated  relative  intensity  of 
the  corresponding  transition  for  the  two  (or  three)  channel 
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FIGURE  6.16 


Vertical  lines  show  calculated  relative  intensities. 


3166 


6000 


3000 


case,  listed  in  Table  four.  Although  it  has  been  establish 
ed  from  reflectance  studies'*''*' ^  that  the  strong  band 

seen  in  transmission  has  a  greatly  broadened  contour,  a 
crude  estimate  of  the  areas  under  the  C&T  bands  at  1450 
and  2044  cm  ^  gives  a  ratio  1:0.34  (±20%),  which  is  in 
good  agreement  with  the  intensity  ratio  calculated, 

1:0.315  (Table  four). 

The  calculation  of  the  positions  and  intensities  of 
these  lines  thus  provides  a  strong  confirmation  of  their 
assignment  to  the  series  The  peak  at  1225  cm  ^ 

has  been  assigned  to  the  bending  mode  V£/  and  the  remain¬ 
ing  lines,  at  1833  cm  ^  and  2426  cm  ■*" ,  are  evidently  the 
combinations  \>  ^  +  j  =  1/2. 

In  the  lower  curve  of  figure  6.16  (3000-6200  cm  "*■)  , 
the  vertical  lines  correspond  to  relative  intensities 
calculated  using  one  channel.  It  is  clear  from  the  cal¬ 
culated  line  positions  of  Table  three,  and  relative 
intensities  of  Table  four,  that  the  lines  at  5100,  5650, 
and  6185  cm"'*'  are  combinations  of  with  the  second 
overtone  of  v^,  that  is,  the  series:  1 3  v  3  *  +  j'v^' 

(j  =  0,1,2).  Experimentally,  it  can  be  seen  that  the  inten 
sities  of  '3v3'  and  '  3v3'  +  ,v1*  are  not  comparable  in 
magnitude.  Table  four,  however,  shows  that  if  the  Fermi 
resonance  effect,  between  |0,3>  and  |  7 , 1> ,  were  present, 
these  intensities  would  be  of  comparable  magnitude,  and 
it  is  concluded  that  the  close  degeneracy  of  the  states 


■ 
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is  an  artifact  of  the  potential  model. 

The  relative  intensity  of  the  lines  at  3166  cm  1  and 

5100  cm  ^  offers  strong  support  for  the  labelling  +  3'v^' 

and  '3v3'  respectively,  in  agreement  with  the  assignment 

of  C&T.  There  has  been  some  disagreement  about  this 

102 

assignment.  Newman  and  Badger  have  assigned  lines  at 

4240  and  4890  cm  to  '  3v^ '  and  '3v3'  +  'v-^'  respectively, 

and  subsequently  5085  and  5590  cm  \  to  +  2v^  +  2v^ 

and  +  2V2  +  3v^.  On  the  basis  of  a  consideration  of 

9  5  2  0 

the  deuterium  spectrum  (KDF2)  of  Ketelaar,  Ibers  has 
argued  in  favour  of  the  C&T  assignment.  This  calculation 
strongly  supports  his  conclusion. 

A  line  at  3000  cm  1  has  been  assigned  to  V2  +  *3v^' 
by  Dawson  et  al . and  the  line  at  about  3400  cm  1  is 
probably  due  to  \>2  +  '  4v^'. 

The  absorption  region  3600-5000  cm  1  is  evidently  due 
primarily  to  overtones  and  combinations  of  the  bending 
modes  (C&T  discuss  this  in  detail) ,  although  the  combina¬ 
tion  +  ' 4v^'  will  clearly  contribute  some  intensity 
in  the  region  of  3700  cm  .  Higher  combinations, 
v3  +  3’vi’  =  5,6,..),  will  not  contribute  significantly 

in  this  region,  except  perhaps  through  Fermi  resonance 

with  the  bending  overtones  and  combinations. 

2  4 

The  comparison  by  Alml^f,  of  the  lines  he  calculated 
at  3916  and  4750  cm”1,  with  experimental  lines  in  the 
region  3800-4800  cm-1  is  invalid,  for  the  following  reasons. 


■ 
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Since  his  potential  surface  (being  only  a  function  of  the 
linear  coordinates  and  C^)  can  only  estimate  stretch¬ 
ing  frequencies,  his  assignment  implies  the  experimental 
lines  are  'stretching'  in  origin,  whereas  it  has  been 
shown  here  that  the  relative  intensities  for  these  stretch¬ 
ing  transitions  are  much  too  low  to  account  for  the 
intensity  observed  in  this  region.  Alml^f  did  not  calculate 
relative  intensities. 

C.  Conclusions. 

1.  Stepanov  Approximation. 

In  the  absence  of  near-degeneracy  between  adiabatic 

protonic  states  (e.g.,  and  2\>2  treated  by  Witkowski 
15 

and  Wojcik  ) ,  the  Stepanov  approximation  has  been  shown 
to  give  acceptable  accuracy  (~1%)  for  the  prediction  of 
vibrational  eigenvalues.  The  inclusion  of  non-adiabatic 
coupling  between  adjacent  adiabatic  protonic  states  (2 
channel  coupling  for  most  IR  transitions)  yields  essentially 
exact  solutions. 

It  does  not  appear  to  be  necessary  to  invoke  the 
mechanism  of  "electrical  anharmonicity "  explicitly  in 
order  to  explain  the  abnormally  high  relative  intensities 
of  combination  bands  and  overtones  which  together  produce 
the  unusually  broad ,  structured  "v^  band"  so  characteristic 
of  the  IR  spectrum  of  an  H-bonded  species.  Relative 
intensities  calculated  within  the  Stepanov  approximation, 
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and  with  a  zero  electrical  anharmonicity  parameter,  are 
of  the  same  order  of  magnitude  as  those  observed  experi¬ 
mentally.  Inclusion  of  adjacent  state  non-adiabatic 
coupling  may  lead  to  about  a  10%  increase  in  combination 
band  intensity,  again  in  the  absence  of  degeneracy.  The 
combination  bands  have  thus  been  shown  to  occur  as  a 
vibrational  analogue  of  the  Franck-Condon  principle. 

For  the  case  of  isolated  "Fermi  resonance",  | 7 , 1 > 
and  |  0 , 3 >  in  the  model  problem,  it  is  clearly  necessary 
to  go  beyond  the  Stepanov  approximation.  A  perturbation 
treatment  in  the  two-dimensional  subspace  could  not  have 
predicted  the  energy  level  reversal  of  order  which  is 
obtained  by  exact  solution  with  the  non-adiabatic  couplings, 
and  of  course  the  Stepanov  eigenvalues  and  wavef unctions 
are  seriously  in  error  for  this  case. 

These  conclusions  are  not  restricted  to  the  particular 
model  problem  (HF2~)  investigated,  because  this  model 
represents  the  strongest  case  of  non-adiabatic  coupling 
likely  to  occur  in  any  hydrogen  bond:  the  potential 
V^(z;R)  has  strong  anharmonic  interaction  of  the  two 
stretching  coordinates,  leading  to  strong  deformation  of 
the  protonic  states  with  R;  the  reduced  mass  ratio  (m/y) , 
and  adiabatically  separated  frequency  ratio  /  are 

unusually  large. 

It  seems  likely,  therefore,  that  the  Stepanov 
approximation  provides  a  useful  basis  for  the  interpretation 
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of  the  vibrational  spectra  of  most  H-bonded  systems  (where 
no  degeneracy-mediated  couplings  occur) .  Even  for  the 
cases  of  isolated  "Fermi  resonance",  where  non-adiabatic 
couplings  must  be  included,  the  adiabatic  representation 
still  provides  a  significant  advantage  over  the  "matrix 
expansion"  methods  (Chapter  one,  section  B2)  in  terms  of 
the  clarity  of  interpretation  (the  "dominant  channel" 
concept) . 

2.  Vibrational  Spectrum  of  the  Bifluoride  Ion. 

Absolute  values  for  and  v^,  for  an  isolated  HF 

ion,  are  not  available.  Experimental  measurements  of 

vary  from  1425  to  1600  ^3  cm  ^  depending  upon  the 

cation  (K+,  Na+  or  NH^+  generally),  and  the  environment 

103 

(Ketelaar  et  al .  have  measured  -  1599,  1570,  1527, 

1478  cm  1  for  KHF2  in  NaBr,  KCl,  KBr,  KI  respectively; 

Salthouse  and  Waddington^^  have  ~  1560  cm  ^  for  KHF2 

in  KCl) .  Also  v ^  has  been  observed  at  600  and  610 

cm’’*'  97,106  110  £or  KHF^  (with  a  correlation  field  splitting 

due  to  the  crystal  site  symmetry) ,  and  at  -630  cm  ^  for 

105 

NaHF2  (no  splitting) .  Azman  and  Ocvirk  have  also 
reported  v-^  -  675  cm  ^  for  NaHF2* 

The  values  calculated  in  this  work  (v^  ~  1520  cm  ^ ; 
v-^  -  677  cm  ^)  strongly  suggest  that  the  potential  function 
used  is  at  least  a  close  model  of  the  physical  reality. 

Keeping  in  mind  that  the  weak  dependency  of  the 
stretching  eigenvalues  on  the  state  of  excitation  of  the 
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bending  mode  has  been  ignored,  the  interpretation  of  the 
IR  spectrum  of  crystalline  KHF2  essentially  confirms  that 
of  Cote  and  Thompson, 

Though  the  bending  modes  have  not  been  explicitly 
calculated  here,  it  is  possible  to  discuss  them  on  the 
basis  of  the  separation  given  in  Chapter  two,  section  Cl 
(equation  2.15) . 

Unfortunately,  the  experimental  situation  with 
respect  to  is  rather  unclear.  Cote  and  Thompson'^'*' 
appear  to  have  observed  a  ~50  cm  ^  splitting  (static  field 
splitting)  in  KHF2  between  bending  frequencies  perpen¬ 
dicular  (1225  cm  ^)  and  parallel  (1274  cm  ^)  to  the 

tetragonal  axis,  but  the  polarized  single  crystal  trans- 

102 

mission  spectra  of  Newman  and  Badger  show  the  splitting 
to  be  only  10  cm  ^ .  Analyses  of  the  reflectance 
spectrum^'*"'*' '  of  KHF2  have  given  a  splitting  of  4-7  cm  ^ , 
with  V2  ~  1227  and  1235  cm  \  In  the  sodium  salt  there 
is  no  splitting,  and  the  frequency  is  only  slightly 
smaller  :  1210  to  1223  cm  \  while  Ault^^  finds 

1217  cm  1  for  HF2-  in  an  argon  matrix  (CsF/HF) .  Cooke 

112  -1 
et  al .  have  reported  a  splitting  of  25  cm  for 

isotopically  isolated  ^2“  in  KDF2  (and  17  cm  ^  for  DF2” 

in  KHF2) /  and  have  interpreted  this  as  a  measure  of  the 

true  local  anisotropy  of  the  bending  force  field  in  the 

lattice  of  the  potassium  salt. 
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As  mentioned  in  Chapter  two,  section  Cl,  however, 
the  model  involving  cylindrical  symmetry  is  easily  modified 
to  allow  for  non-degenerate  bending  modes.  The  essential 
feature  is  the  R-dependence  of  the  force  constants  for 
these  motions.  This  leads,  via  the  analogous  Franck- 
Condon  principle  to  the  series  of  combination  bands 
v2  +  which  is  clearly  present  in  the  spectrum  of 

Cote  and  Thompson'^'*"  (Figure  6.16).  A  monotonic  R-dependence 
of  the  force  constant  would  mainly  displace  the  minimum  of 
the  potential  curve  for  the  first  excited  state  of  the 
bending  motion  (relative  to  the  ground  state  curve) , 
rather  than  strongly  modifying  the  "force  constant  for 
R  motion"  in  that  curve. 

It  is  probable  that  the  moderately  intense  features 

in  the  regions  3660-3830  cm  \  and  4200-4500  cm  \  observed 

by  Cote  and  Thompson"*"^  and  Newman  and  Badger,  are  not 

due  to  the  bands  "3V2"  and  "3V2  +  V-,"  as  C&T  have  supposed. 

The  dipole  moment  integral  in  p  space  (analogous  to 

equation  (5.33),  for  I  (R)  in  z  space)  for  the  second 

KJ 

overtone  transition  t  has  appreciable  value  at  any 

R  point  only  if  the  potential  function  for  bending  motion 
is  strongly  anharmonic  in  p.  Since  this  is  unlikely,  it 
would  be  expected  that  the  intensity  of  " 3V2"  should  be 
very  much  less  than  the  intensity  of  "3^"  (for  which  the 
potential  function  in  z  space  is  highly  anharmonic) ;  yet 
the  band  at  3660-3830  is  in  fact  stronger  than  "Sv^" . 
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An  alternative  assignment  is  possible.  Cooke, 

112 

Pastorek,  Carlson  and  Decius  have  pointed  out  that  m 
KHF2  (but  not  in  NaHF^)  intermode  coupling  between  \>2  and 
may  occur  (they  have  components  of  the  same  unit  cell 
symmetry) .  It  is  possible  that  the  potential  surface 
(for  bending  in  one  molecule,  and  simultaneous  protonic 
stretching  in  its  neighbour)  could  have  significant 
anharmonic  interaction  terms,  leading  to  the  appearance 
of  the  bands  "2v 2  +  and  +  ^2"  with  relatively 

high  intensity. 

It  is  therefore  proposed  that  the  band  at  3600-3830 
cm  ■*"  should  be  associated  with  the  transition  "2V2  +  v^"/ 
and  that  at  4200-4500  cm  ^  with  +  v2" •  If  such  an 

intermolecular  interpretation  is  correct,  then  (i)  the 
bands  should  not  appear  in  the  spectrum  of  NaHF2  (where 
the  unit  cell  contains  only  one  HF2-  ion) ,  and  (ii)  isotopic- 
ally  isolated  spectra  of  HF2-  in  KDF2  (and  DF2-  in  KHF2 ) 
should  give  different  locations  for  these  bands  (compared 
to  pure  salts) . 
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